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ABSTRACT 

Laser  radar  analyses  up  to  now  generally  fall  into  two  categories.  Single 
pixel  analyses  exist  which  cover  receiver  design,  target  statistics, 
atmospheric  effects,  and  the  resulting  statistical  performance 
characterizations  like  probabilities  of  detection  (FL)  and  false  alarm 
(Pp).  Multipixel  analyses  exist  which  use  various  aa  hoc  target 
recognition,  identification,  and  feature  extraction  schemes  based  on 
intuitive  insight  or  analogies  to  human  perception.  The  limitation  of 
these  multipixel  processors  is  it  is  difficult  to  quantitatively  predict 
processor  performance  as  a  function  of  various  system  parameters  and  the 
resolving  power  in  the  measurement  dimensions  like  range,  Doppler  shift, 
and  angle.  The  single  pixel  performance  equations  can  give  quantitative 
answers  to  these  questions,  but  only  for  single  pixel  measurements.  This 
thesis  bridges  these  two  approaches  by  proposing  physically  realistic 
target,  background,  and  radar  models  which  allow  us  to  incorporate  the 
statistics  for  the  single  pixel  measurements  into  multipixel  probability 
density  functions  and  derive  quasi-optimal  generalized  likelihood  ratio 
processors  from  the  densities. 

^  - - 

The  thesis  proposes  simple  but  realistic  models,  then  derives  the 
multipixel  density  functions,  the  quasi-optimal  multipixel  processors,  and 
Fd-Pf  curves  for  the  processors  as  a  function  of  the  radar  parameters  and 
measurement  dimension  resolving  powers.  The  analyses  are  for  a  ranging 
radar,  but  duality  arguments  are  presented  which  demonstrate  how  to  apply 
the  results  directly  to  Doppler  and  range  and  Doppler  radars.  Some 
analyses  are  presented  which  show  how  the  radar's  resolution  in  range  and 
angle  measurement  dimensions  depend  on  each  other  through  the  radar 
equation  for  the  carrier-to-noise  ratio,  how  they  interact  to  affect 
performance,  and  how  it  is  possible  to  balance  performance  loss  in  one 
measurement  dimension  by  improved  resolution  in  another. 
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I.  INTRODUCTION 


A.  Laser  Radar  Overview 

The  development  of  lasers  with  sufficient  stability,  spectral  purity,  and 
power  has  allowed  system  designers  to  translate  much  of  microwave  radar 
design  directly  to  the  optical  wavelength  region  [1].  Because  their 
wavelength  is  so  much  shorter,  laser  radars  offer  spatial,  angular,  range, 
and  Doppler  shift  resolutions  far  superior  to  those  available  with 
microwave  systems.  Although  laser  systems  can  take  advantage  of  their 
increased  resolution  in  any  or  all  of  these  dimensions,  the  increased 
spatial  and  angular  resolutions  have  particularly  affected  the  way 
designers  present  data  to  the  user.  With  beam  diameters  of  a  few 
centimeters  at  the  radar  optics  exit  aperture  and  beam  angles  measured  in 
microradians,  it  is  possible  to  scan  the  beam  over  the  target  and 
construct  an  image  of  the  target,  similar  to  a  television  (TV)  image,  from 
the  radar  returns.  The  color  and  intensity  of  the  displayed  image  can 
represent  reflected  intensity,  Doppler  shift,  range,  or  some  combination 
of  these  signal  dimensions  [2]. 

With  the  shorter  wavelength  of  laser  systems,  there  comes  a  disadvantage 
as  well.  Because  the  wavelength  is  so  short,  most  target  surfaces  are 
rough  at  dimensions  on  the  order  of  the  radar  wavelength.  So,  they 
reflect  the  incident  beam  with  an  essentially  random  phase  shift  as  a 
function  of  position  on  the  target,  producing  a  phenomenon  known  as  laser 
speckle  [3].  Laser  speckle  appears  as  a  graininess  in  the  image  of  the 
target  due  to  self  interference  among  the  reflected  wave's  angular 
components.  Since  it  is  a  self  interference,  it  is  independent  of  signal 
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strength.  It  appears  even  at  very  high  carrier  (reflected  radar  energy) 
to  noise  ratios.  Any  processing  scheme  should  take  into  account  both 
detection  noise  and  speckle  noise,  both  of  which  degrade  image  quality. 

B.  Image  and  Signal  Processing 

Laser  radars  (and  some  microwave  radars  as  well)  have  sufficient  spatial 
resolution  that  they  can  build  images  of  targets  by  making  spatially 
distributed  measurements  of  the  target  with  a  beam  of  diameter  much  less 
than  the  target  dimensions.  Each  spatially  distinct  element  of  the  target 
is  called  a  picture  element,  or  pixel.  Up  to  now,  analyses  for  these  high 
resolution  radars  have  fallen  into  two  general  classes:  single  pixel 
statistical  signal  processing  and  multipixel  image  processing. 

The  single  pixel  processing  work  follows  the  traditional  radar  and 
communication  system  approaches  [4],  This  approach  models  the  target  and 
the  radar,  derives  statistically  optimal  processors  based  on  the  models, 
and  produces  equations  for  the  performance  of  the  processors  derived  from 
the  statistics.  The  advantages  of  this  approach  are  twofold.  First,  the 
performance  equations  are  in  terms  of  the  various  laser  system  design 
parameters.  This  makes  it  easy  to  predict  how  performance  will  change  as 
we  vary  the  system  parameters.  Second,  since  the  processors  are 
statistically  optimal,  their  performance  measures  must  be  bounds  on  the 
performance  achievable  with  any  other  proposed  processor.  The 
disadvantage  of  this  approach  is  it  is  a  single  pixel  analysis.  The 
processors  don’t  use  the  information  available  to  them  from  the  spatially 
distributed  measurements  of  several  pixels. 


The  image  processing  approach,  on  the  other  hand,  often  draws  on  pattern 
recognition  and  image  analysis  theory  from  research  in  areas  like  remote 
sensing,  machine  vision,  and  artificial  intelligence  [ 5 ] ,  [6].  This 
approach  proposes  intuitive  algorithms  based  on  these  research  fields, 
and  the  resulting  ad  hoc  processors  make  little  use  of  the  statistics  of 
the  signals  generating  the  images.  The  advantage  of  this  approach  is  many 
of  the  processors  seem  to  work  quite  well.  The  disadvantage  is  it  is 
incapable  of  producing  performance  equations  which  allow  system  designers 
to  predict  how  changes  in  the  radar  parameters  will  affect  performance. 

This  thesis  presents  for  the  first  time  a  method  to  perform  statistical 
signal  processing  based  on  target  and  laser  models  like  the  single  pixel 
approach  just  outlined,  but  which  takes  advantage  of  the  multipixel  nature 
of  the  targets.  The  result  will  be  equations  for  the  radar's  performance 
against  spatially  distributed  or  multipixel  targets  in  terms  of  the 
radar's  design  parameters.  These  equations  are  important  because  they 
give  designers  a  method  to  perform  trades  among  system  parameters  and 
predict  the  effect  on  system  performance.  With  these  equations  it  is 
possible,  for  instance,  to  trade  intensity  resolution  with  range 
resolution  or  spatial  resolution  with  Doppler  resolution  and  predict  the 
effects  on  receiver  performance.  It  is  not  possible  to  do  this  directly 
with  the  ad-hoc  models.  These  equations  also  act  as  benchmarks  against 
which  we  can  compare  the  performance  of  the  ad  hoc  image  processing 
algorithms  researchers  will  undoubtedly  continue  to  pursue. 

C.  Problem  Statement 

We  will  study  the  classical  radar  detection  problem  [4],  In  our  detection 
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problem  we  will  assume  we  are  given  a  set  of  directions  comprising  an 

angular  uncertainty  region,  ftu,  and  a  corresponding  swath  of  ranges 

representing  a  range  uncertainty  region,  Lu,  wherein  a  target  may  appear, 

as  in  Figure  1.1.  There  may  or  may  not  be  a  target  in  this  uncertainty 

volume  ft  L  and,  if  there  is,  we  don’t  know  where  it  falls  within  the 
u  u 

volume.  Our  task  is  to  make  measurements  over  the  uncertainty  volume  with 
the  laser  radar  and,  based  on  the  measurements,  decide  between  two 
hypotheses:  H^,  target  present,  and  Hq,  target  absent.  Although  we  will 
try  to  treat  thoroughly  general  processors,  the  goal  here  is  to  derive 
tools  to  predict  the  performance  of  real  systems.  If  necessary,  we  will 
make  simplifying  assumptions  to  get  results  which  are  useful  for  existing 
systems  at  the  expense  of  complete  generality  in  the  analysis. 


D.  Thesis  Organization 

Chapter  two  of  this  thesis  develops  background  information  on  laser  radar 
systems  and  target  statistical  models.  It  derives  a  statistical  model  for 
the  signal  in  the  intermediate  frequency  (IF)  system  of  the  radar  and 
imposes  some  IF  pre-processing  used  in  real  laser  radar  systems.  This 
chapter  also  takes  the  crucial  step  of  presenting  a  model  for  the  major 
characteristics  of  the  target  and  its  environment  in  a  fashion  which  makes 
it  possible  to  incorporate  these  characteristics  into  the  statistical 
signal  model. 


Chapter  three  derives  a  joint  probability  density  function  for  two  of  the 
most  common  radar  measurements:  range  and  intensity.  Chapter  three 
closes  by  considering  the  range-Doppler  and  range-range  and  Doppler 
dualities  which  make  it  possible  to  do  the  Doppler  and  range  and  Doppler 
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system  analyses  by  direct  analogy  with  the  range  system  analysis. 


Chapter  four  introduces  the  binary  detection  receiver  and  the  generalized 
likelihood  ratio  and  derives  optimal  and  near  optimal  processors  for 
several  measurement  models.  Chapter  five  takes  two  of  the  most  important 
processors  from  chapter  four  and  derives  equations  to  characterize  their 
performance. 

Chapter  six  presents  numerical  and  graphical  results  for  the  equations  in 
chapter  five.  These  results  demonstrate  how  the  radar  performance,  in 
terms  of  probability  of  detection,  P^,  and  probability  of  false  alarm,  Pp, 
varies  as  a  function  of  the  laser  design  parameters,  the  target 
characteristics,  and  the  type  of  radar  used.  Specifically,  we  will  assess 
how  performance  varies  as  a  function  of  radar  resolution  capabilities  for 
several  types  of  radars  with  laser  sources  of  either  constant  average  or 
peak  power.  The  ultimate  output  of  this  analysis  is  the  ability  to 
predict  for  the  first  time  what  sets  of  design  parameters  will  allow  the 
radar  to  meet  a  given  performance  level  against  a  multipixel  target. 

The  final  chapter  summarizes  the  thesis,  draws  some  conclusions,  and 
outlines  some  areas  where  further  research  may  be  profitable. 


II.  USER  RADAR  BACKGROUND 


A.  System  Configuration 

Figure  2.1  is  the  simplified  block  diagram  of  a  coherent  laser  radar. 

This  thesis  considers  only  monostatic  radars  —  the  transmitter's  exit 
optics  and  the  receiver's  entrance  optics  are  the  same.  This  could  be  a 
ranging  or  Doppler  radar  (or  both).  The  principal  differences  between 
these  systems  are  the  transmitted  waveforms  and  the  type  of  intermediate 
frequency  (IF)  processing.  For  a  ranging  system,  the  waveform  is  usually 
as  short  as  possible  to  give  good  range  resolution.  For  a  Doppler  system, 
the  waveform  may  be  continuous  wave  (CW)  or  a  long  pulse.  It  is  possible 
to  time-share  both  measurements  with  an  appropriate  choice  of  waveform 
[7].  The  radar  performs  a  raster  scan  of  the  transmitter  beam  as  shown  in 
figure  2.2.  The  scan  rates  are  such  that  there  is  minimal  beam  overlap 
between  adjacent  pixels.  The  radar  may  be  on  a  truck,  an  aircraft,  or 
elsewhere.  For  the  air-to-ground  imaging  situation,  the  radar  is  above 
and  looking  down  on  the  target.  This  imposes  a  specific  laser-target 
geometry,  which  we  will  introduce  and  exploit  in  section  E  of  this 
chapter. 

This  thesis  uses  the  laser  radar  systems  in  existence  at  MIT  Lincoln 
Laboratory  as  examples  of  real  systems  when  desiring  to  make  theoretical 
models  simpler  or  more  realistic,  or  to  use  specific  numerical  values  for 
design  parameters.  The  Lincoln  Laboratory  systems  use  CC^  lasers 
operating  at  a  wavelength  of  10,6  ym.  These  systems  operate  at 
relatively  short  ranges  (on  the  order  of  a  few  kilometers)  with  modest 
aperture  sizes  (on  the  order  of  10  to  20  centimeters)  [8],  [9].  The 
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Figure  2.2 

Radar  Beam  Raster  Scan  in  Angle  Space 
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aperture  size,  laser  power  level  (up  to  a  few  watts  [8],  [9]),  and 
detector  sensitivity  combine  to  limit  the  radar's  range  under  any  given 
set  of  atmospheric  conditions.  The  aperture  size  is  not  larger  because 
these  systems  are  designed  to  investigate  the  design  of  primarily  compact 
systems  for  air-to-ground  applications.  These  systems  also  use  heterodyne 
detection  rather  than  direct  detection.  The  advantages  of  a  heterodyne 
system  are  that  it  increases  the  sensitivity  of  the  detector  because  of 
the  well  known  detection  gain  of  a  heterodyne  system  [10]  and,  since  the 
local  oscillator  is  coherent  with  the  received  optical  signal,  it 
preserves  the  relative  phase  variations  of  the  received  signal  so  the 
receiver  can  measure  Doppler  shifts  in  the  optical  carrier  [10],  Since  a 
direct  detection  receiver  is  not  sensitive  to  the  phase  of  the  optical 
carrier,  it  cannot  measure  Doppler  shifts.  The  disadvantage  of  a 
heterodyne  system  is  the  greater  prominence  of  the  speckle  effect.  In  a 
direct  detection  system,  it  is  possible  to  reduce  the  degree  of  speckle 
induced  intensity  variations  by  using  a  larger  diameter  detector  to 
average  over  several  speckle  lobes.  In  a  heterodyne  system,  the 
diffraction  limited  field-of-view  makes  it  impossible  to  reduce  speckle  so 
easily  [11]. 

To  analyze  this  configuration,  we  will  first  model  the  target,  then  use 
Fourier  optics  to  propagate  the  laser  pulse  to  and  from  the  target,  and 
use  tools  from  statistical  communication  theory  to  determine  the 
statistics  of  the  signal  just  before  IF  processing  for  a  single  pixel. 

Then  we  will  discuss  the  pre-processors  often  used  in  laser  radar  systems 
and  their  advantages  and,  finally,  present  a  model  for  the  entire 
multipixel  target  and  its  surroundings. 
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B.  Target  Models 

Real  targets  are  themselves  very  complicated  objects.  However,  this 
thesis'  purpose  is  not  to  investigate  the  subtle  differences  between  two 
trucks  from  the  same  assembly  line,  or  two  different  trucks,  or  even  a 
truck  and  a  building.  Accurate  statistical  modeling  of  real  targets  is 
the  subject  of  other,  ongoing  research  [12].  In  this  thesis,  all  targets 
will  follow  the  purely  diffuse  speckle  target  model  [3].  This  model  is 
very  accurate  for  certain  prepared  test  targets  [13],  [14]  and  quite 
accurate  for  at  least  the  gross  features  of  real  targets.  In  the  most 
general  form  of  this  model,  the  target  has  a  complex  reflection 
coefficient  IX  p,  z,  t  )  at  every  point  in  space  and  time  (  p,  z,  t  ), 
where  p  is  the  transverse  position  vector  (  x,  y  )  and  T(  *  )  is  a  four 
dimensional  (space  and  time)  complex  Gaussian  random  process  with 
independent,  identically  distributed  real  and  imaginary  parts  [3].  Since 
T(  *  )  is  Gaussian,  a  second  moment  description  completely  characterizes 
its  statistics.  The  moments  are: 


E[  T(  p,  z,  t  )  ]  =  0, 

E[  T(  P,  z,  t  )  T(  p',  z\  tr  )  ]  =  0, 

E[  T(  p,  z,  t  )  T*(  p\  z\  t’  )  ]  = 


(  2.1  ) 


X  Tg(  p,  z,  t  -  t'  )  <5(  p  - 


o')  S(z-z'). 


The  spatial  impulse  functions  indicate  the  random  structure  of  the 
target's  surface  decorrelates  much  faster  than  a  distance  of  one 
wavelength  X  at  the  laser  frequency,  Vq  =  c  /  X  .  The  spatial 
dependences  in  T  (  *  )  indicate  the  target  need  not  be  uniform.  The  time 
dependence  indicates  the  target  may  be  moving,  but  this  model  assumes  the 
motion  is  a  stationary  random  process,  hence  the  dependence  on  only  the 


-  20  - 


time  difference  t  -  t'.  This  will  be  the  model  (often  simplified)  for  all 

targets,  background,  clutter,  or  anything  entering  the  radar’s  field  of 

view  (fov),  including  the  no  target  case  where  T  (  p  ,  z,  t  )  =  0. 

s 


C.  IF  Signal  Derivation 

In  terms  of  the  complex  envelope  of  the  optical  signals  at  the  optical 
frequency,  v  ,  the  IF  signal  out  of  the  detector  is  proportional  to  the 
product  of  the  complex  field  envelope  of  the  reflected  pulse  and  the  local 
oscillator  (L0)  integrated  over  the  detector  surface.  To  get  the 
reflected  pulse's  complex  envelope,  one  must  propagate  the  transmitted 
field  to  the  target,  multiply  by  the  target  reflection  coefficient  T(  *  ) 
and  propagate  back  to  the  receiver.  In  this  thesis  we  will  model 
atmospheric  optical  propagation  as  free-space  diffraction  modified  only  by 
extinction  loss,  which  is  a  reasonable  approach  for  compact  CO^  laser 
radars  [24],  We  shall  assume  the  L0  field  pattern  is  matched  to  the 
transmitted  field  pattern,  as  is  usually  the  case  in  such  radars  [14],  In 
practice,  it  is  easier  to  use  the  antenna  theorem  [15]  and  back  propagate 
the  L0  field  to  the  target  and  perform  the  integration  over  the  target 
surface  rather  than  the  detector  surface.  Define  £ (  p,  z,  t  )  as  the 
normalized  transmitted  beam  pattern  at  target  range  z.  Here  the  time 
dependence  t  is  due  only  to  the  beam  scanning  motion.  We  define  another 
term,  js(  t  ),  to  contain  the  time  dependence  from  the  transmitter 
modulation.  Now  we  can  express  the  transmitted  beam  and  back  propagated 
L0  beam  at  the  target  in  terms  of  £(  *  )  and  sji  t  ).  Further  define 
and  A^,  as  the  optics  aperture  and  target  surface,  respectively,  and  P^,  as 
the  transmitted  power.  Then  the  complex  envelope  of  the  IF  signal  at 
intermediate  frequency  (radians/second)  is  t  )  [14]: 


-  21  - 

i(  t  )  =  /  dz  s(  t  -  2z/c 

•T(o,z,  t-z/c)  g 

where  n{  t  ),  which  represents  the  LO  induced  shot  noise,  is  a  stationary, 
zero-mean,  circulo-complex  Gaussian  random  process  independent  of  X(  *  ) 
with  covariance  function  E[  ji(  t  )  ri*(  t’)]  =  (hv  /n)  6  (  t  —  t '  ) 

[10],  Here  h  vq  is  the  optical  photon  energy  and  n  is  the  detector 
quantum  efficiency.  The  first  term  in  equation  2.2  is  basically  a  delay- 
integral/heterodyne-overlap  integral  specification  of  the  target-return 
component  of  j[(  t  ).  The  delay  integral  over  z  encompasses  targets  at  all 
ranges  while  the  overlap  integral  over  A,p  is  the  heterodyne  integral  for  a 
particular  range  z.  The  equation  includes  all  appropriate  lag  times.  For 
the  IF  target  return  at  time  t  from  range  z,  the  transmitted  waveform  was 
_s(  t  —  2z/c  ).  The  first  C  (  *  )  term  is  the  transmitted  beam  shape  that 
contributes  to  the  target  return  at  t  from  range  z.  The  second  _§_(  *  ) 
term  is  the  LO  beam  shape  relevant  to  the  target-return  at  t  from  range  z; 
it  has  no  time  lag.  The  target  term  T.(  *  )  has  a  lag  of  only  z/c  seconds, 
the  one-way  time  of  flight  from  the  target  to  the  radar. 

Now  let  us  make  the  assumption  g  (  p  ,  z,  t-2z/c  )  -  p"  ,  z,  t  )  for 

all  z  such  that  I  £  (  *  )|  makes  an  appreciable  contribution  to  the 
integral.  This  is  the  assumption  of  negligible  radar  lag  angle.  It  will 
make  more  physical  sense  after  defining  the  transmitted  field  pattern  in 
terms  of  a  radar  scanning  motion  in  the  next  paragraph.  For  now,  the  IF 
waveform  is: 


)  I  dp  _£_(  P  ,  z,  t  -  2z/c  ) 
(p,z,t)  +n(t),  (  2.2) 
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|  .  00  ^ 

jr(  t  )  =  J  dz  s(  t  -  2z/c  )  J  dp  T(  p  ,  z,  t  -  z/c  ) 

0  Aip 


2 

_L_  (  P  ,  Z,  t  -  z/c  )  +  n(  t  ). 


(  2.3  ) 


Although  the  transmitted  beam  pattern  has  been  arbitrary  so  far,  it  is 
usually  a  collimated  Gaussian  beam  with  some  pointing  angle  d>  (  t  )  with 
respect  to  a  reference  direction  (see  figure  2.4).  Thus,  let  us  take 
H(  "P  ,  t  ),  defined  to  be  the  normalized  transmitted  beam  pattern  at  the 
radar,  as  follows: 

2 

F(  p-  ,  t  )  =  (  )1^2  exp[  -  ^  I  I  +  jk  “p  •?(  t  )  ]  (  2.4) 

^  aa  al 


where 


k  =  2  tt  /  \  ~  wavenumber 
-2 

d  =  e  beam  intensity  diameter  at  exit  pupil. 


Generally,  <j>  (  t  )  is  a  linear  scanning  function,  1j>(  t  )  =  4>(  tn  )  + 

4  •  (  t-t  ),  for  t  <  t  <  t  +  T.  ,  where  t  and  "$(  t  )  are  the  starting 
time  and  direction,  respectively,  for  the  n-th  scan  line,  $  is  the  scan 
rate,  and  is  the  time  required  to  scan  one  line.  This  gives  a  raster 
scan  arrangement  like  that  in  figure  2.2.  (Note  that  because  of  flyback 
time,  tn+^  >  tn  +  T^  generally  prevails). 


D.  IF  Pre-processing 

Thus  far  we  have  made  as  few  assumptions  as  possible  in  order  to  make  the 
equation  for  _j_(  t  )  as  general  as  possible.  These  derivations  give  the 
statistics  of  the  IF  waveform  directly  out  of  the  detector,  but  real 
processing  systems  (i.e.,  detection,  identification,  etc.)  generally  don't 


-  J  J**  j*  j*  -f-  ^  j\  *  *  r  *•  ' 
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deal  with  this  IF  waveform  directly.  Instead,  some  form  of  pre-processing 
usually  exists  between  the  IF  and  the  signal  processing  [12],  [13].  One 
reason  for  this  intermediate  step  is  the  bandwidth  disparity  between  the  IF 
waveform  and  the  bandwidth  of  recording  and  processing  electronics. 

Typical  research  systems  record  the  radar  return  data  on  tape  for 
subsequent  off  line  analysis  and  testing  [12],  [13].  It  is  difficult  to 
record  the  data  coming  directly  off  of  the  detector  because  the  IF 
frequency  is  in  the  tens  of  MHz  and  has  a  bandwidth  of  several  MHz. 

Reducing  this  bandwidth  is  a  prime  reason  why  these  systems  employ  pre¬ 
processing  hardware  prior  to  data  recording.  Reducing  the  quantity  (and 
hence  the  bandwidth)  of  the  data  also  makes  it  easier  to  process  the  data 
in  real  time,  an  important  consideration  for  practical  systems. 


Scanning  systems,  such  as  the  one  we  have  specified,  usually  break  the 
target  into  pixels  and  build  an  image  of  the  target  area  from  the  pixels. 
The  pre-processing  electronics  act  on  each  individual  pixel  time  element. 
For  a  ranging  radar,  the  usual  pre-processor  examines  the  pixel  time 
interval  and  selects  the  peak  of  the  random  waveform  and  records  its 
intensity  and  time  location  relative  to  the  start  of  the  interval 
(corresponding  to  the  target's  range).  For  a  Doppler  radar  system,  the 
usual  pre-processor  Fourier  transforms  the  received  random  waveform  and 
then  selects  the  peak  in  the  frequency  space  and  records  its  intensity  and 
location  (target  velocity).  These  peak  detection  pre-processors  reduce  the 


entire  random  waveform  over  a  time  interval  to  two  real  random  variables. 
Figure  2.3  is  a  schematic  diagram  of  both  pre-processors. 


If  the  pixel  dwell  times  are  short  compared  to  the  time  orders  on  which  the 
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target  changes,  there  is  little  loss  in  optimality  from  dividing  the 
waveform  into  pixels,  but  the  peak  detection  may  seem  a  bit  ad-hoc. 

Actually,  there  is  theoretical  support  for  this  type  of  pre-processing.  In 
the  case  of  a  ranging  radar  and  a  non-range  spread  target  observed  on  a 
single  pixel,  it  is  easy  to  show  [4]  that  the  maximum  likelihood  (ML) 
estimate  of  the  target  range  corresponds  to  the  peak  of  the  matched  filter 
output  of  the  detector.  So  the  pre-processor  is  essentially  the  ML  range 
estimator,  but  remains  fixed  for  all  types  of  problems. 

One  could  argue  the  virtues  of  working  out  performance  measures  on  truly 
optimal  receivers,  given  one  goal  is  to  examine  the  limits  on  the 
performance  of  laser  radars.  However,  real  radars  are  likely  to  continue 
to  use  some  of  the  mentioned  types  of  pre-processing  for  the  reasons 
mentioned.  Taking  this  pre-processing  into  account  will  give  performance 
measures  which  are  more  realistic  and  closer  to  those  researchers  are 
likely  to  achieve.  This  thesis  will  use  the  matched  filter,  peak  detector 
pre-processor  model  throughout.  The  signal  analysis  for  the  pre-processor 
will  simply  pick  up  where  the  analysis  for  _y_(  t  )  left  off.  None  of  the 
previous  work  goes  to  waste.  Detailed  single  pixel  peak  detector  analyses 
exist  [23].  The  analyses  in  this  reference  include  thresholding  in  the  pre¬ 
processor  and  the  effect  of  the  resulting  dropouts.  (When  the  intensity 
for  a  particular  pixel  is  below  the  threshold,  no  data  is  recorded  for  that 
pixel  and  we  say  a  dropout  has  occurred.)  This  thesis  will  not  use 
thresholding  in  the  pre-processor. 

E.  Frames  and  Subframes 

The  pixel  is  our  basic  unit  of  spatial  measurement  and  we  will  shortly 


derive  the  statistics  for  measurements  taken  on  a  single  pixel.  However, 
laser  radars  have  resolutions  such  that  interesting  targets  often  take  up 
more  than  a  single  pixel.  We  would  like  our  receivers  to  use  this 
spatially  distributed  information  in  an  optimal  or  nearly  optimal  fashion 
to  improve  their  performance.  Since  we  have  arranged  our  scanning  geometry 
(see  figure  2.2)  so  pixels  are  essentially  non-overlapping,  the  pixel 
measurements  will  be  independent  thanks  to  the  statistics  of  the  target 
reflection  coefficient  T(  *  ).  Now  we  need  a  method  to  introduce  the 
spatial  relationships  between  the  pixels  into  the  statistical  models. 

We  will  adopt  the  down-looking  geometry  of  figure  2.4  for  our  analysis  of  a 
single  target  at  an  unknown  range  and  angular  location  within  the  radar's 
field  of  regard  (total  field  of  view).  For  this  geometry,  the  radar  always 
illuminates  something.  It  may  be  the  desired  target  (the  vertical  object 
in  figure  2.4),  or  it  may  be  the  horizontal  background.  In  either  case, 
since  the  radar  illuminates  something,  there  is  a  target,  or  reflector 
element,  return  for  each  pixel.  We  will  sometimes  use  the  term  target 
loosely  to  refe.  to  either  target  or  background  reflector  elements.  We 
will  discuss  figure  2.4  in  more  detail  shortly. 

We  will  assume  we  know  the  target  shape  and  angular  extent  so  we  can  tile 
the  angular  observation  plane  with  M  target  shapes,  each  containing  N 
pixels  as  in  figure  2.5.  We  will  call  one  observation  of  the  entire 
observation  plane  (MN  pixels)  a  frame  and  call  each  target  shape  a  subframe 
(N  pixels).  The  radar  may  stare  in  a  particular  direction  long  enough  to 
take  more  than  one  frame  of  data  (MN  pixels),  but  we  will  consider  only  one 
frame  at  a  time  —  we  won't  try  to  use  information  from  one  frame  to  the 
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next.  We  shall  assume  henceforth  the  target  exactly  aligns  with  one  of  the 
subframes.  This  assumption  will  allow  us  to  search  the  frame  for  the 
target  one  subframe  at  a  time  with  each  subframe  independent  of  all 
others.  Without  this  assumption,  searching  the  frame  for  the  target  would 
involve  searching  approximately  MN  subframes,  each  with  many  pixels  in 
common  with  several  other  adjacent  subframes.  That  is,  we  would  be  sliding 
the  target  shape  as  a  window  around  the  entire  frame  like  the  dotted  target 
outline  in  figure  2.5.  Thus,  the  subframes  would  not  be  statistically 
independent  and  this  would  considerably  complicate  the  ultimate  task  of 
computing  receiver  performance.  This  may  seem  a  rather  weak  justification, 
but  we  must  consider  our  stated  objective  is  to  find  useful  performance 
measures  even  at  the  expense  of  optimality  if  necessary.  This  is  our  model 
of  the  target's  angular  characteristics.  Now  consider  the  target's  range 
characteristics. 


We  will  begin  our  range  model  by  assuming  both  the  target  and  the 
background  are  range  unresolved.  A  range  unresolved  object  is  one  where, 
for  each  pixel,  the  radar  beam  illuminates  only  the  target  or  the 
background  for  each  pixel  (we  ignore  effects  at  the  edges  of  the  targets) 
and  the  target  or  background  has  a  sharply  defined  range  relative  to  the 
radar's  range  resolution.  The  iso-range  lines  on  figure  2.4  are  separated 
by  a  distance  equal  to  the  radar’s  range  resolution.  True  to  the 
assumption,  the  target’s  depth  is  much  less  than  the  radar's  range 
resolution  distance.  Under  this  model,  the  z-dependence  of  T(  p  ,  z,  t  ) 
is  simply  a  Dirac  delta  function  at  the  target  or  background  range  L: 

T(  p  ,  z,  t  )  =  T(  p  ,  L,  t  )  <$(z-L).  Of  course,  this  will  make  the  z 
integral  in  equation  2.3  trivial,  but  we  will  wait  until  the  next  chapter 
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to  actually  perform  this  integration.  Another  type  of  interesting  target 
we  won't  consider  in  this  thesis  is  the  range  spread  target.  Range  spread 
targets  are  those  whose  geometry  is  such  as  to  generate  returns  from  a 
spread  of  distances  in  a  single  pixel.  For  example,  a  flat  plate  viewed  at 
close  to  grazing  incidence  would  produce  a  range  spread  return,  as  do 
aerosols  and  the  edges  of  hard  targets.  The  range  unresolved  assumption  is 
reasonable  for  the  2-D  pulsed  imager  compact  CO^  systems  we  are  considering 
in  this  thesis  [8],  [9]. 


Examining  figure  2.4  again,  notice  the  down  looking  geometry  will  make  the 
background  appear  to  slope  away  from  the  radar  at  a  known  angle  (the 
radar's  pointing  angle).  If  we  know  the  radar's  altitude,  we  can  compute 
the  range  to  the  background  accurately  for  a  smooth  background.  So,  we 
will  assume  we  know  the  range  to  the  background  for  each  pixel  in  the 
frame.  We  will  let  any  target  in  the  frame  be  at  an  unknown  range,  but 
require  it  to  be  vertical.  Most  targets  of  interest  easily  meet  this 
requirement.  For  a  reasonably  smooth  background,  each  horizontal  line  of 
pixels  in  a  subframe  will  be  at  the  same  range  and  adjacent  lines  will  be 
at  different  ranges.  The  target,  however,  will  have  all  its  pixels  in  the 
subframe  at  the  same  range.  Notice  if  the  target  sits  on  the  ground,  as  in 
figure  2.4,  its  range  is  the  same  as  the  range  for  the  nearest  background 
pixel  in  the  subframe.  But  since  we  know  the  range  to  every  background 
pixel,  we  would  know  the  target  range,  too.  However,  the  target  need  not 
actually  sit  on  the  ground;  it  may  be  above  the  ground  as  well.  As  a 
result  we  must  assume  we  don't  know  the  target  range.  However,  since  the 
target  can't  be  behind  the  ground,  we  know  the  background  range  is  always 
more  distant  than  the  target. 


This  completes  our  model  of  the  target's  range  and  angular 
characteristics.  We  have  introduced  a  model  which  is  complex  enough  to  be 
fairly  realistic  while  being  simple  enough  for  analysis.  We  have  made  the 
model  complex  by  allowing  the  target's  angular  position  and  range  be 
unknown,  by  forcing  the  contiguous  nature  of  the  target  into  the  subframe 
model,  and  introducing  a  background  much  like  the  target.  Both  target  and 
background  have  unknown  reflectivities.  The  model  is  simple  by  virtue  of 
the  relatively  small  number  of  free  parameters  in  the  model.  These 
parameters  are  the  target  range  and  reflectivity,  L>t  and  p  t,  the 
background  reflectivity,  p  and  the  single  subframe  containing  the 
target,  mQ.  The  target  size  and  shape  are  known  and  both  background  and 
target  have  uniform  average  reflectivity.  Finally,  all  the  pixels  are 
independent  as  a  result  of  the  speckle  statistics. 

In  the  next  chapter  we  will  derive  probability  density  functions  for  the 
pre-processor  outputs,  which  are  our  measurement  statistics,  starting  from 
the  j[(  t  )  derivations  in  this  chapter.  In  chapter  four  we  will 
incorporate  the  target  frame  models  above  and  the  chapter  three  results 
into  near  optimal  receiver  processor  equations. 


III.  RANGE-INTENSITY  DENSITY  FUNCTION 


In  the  previous  sections  we  were  able  to  describe  the  statistics  of  the 
complex  envelope  of  the  detector's  output  exactly.  It  had  a  Gaussian 
density  function  with  moments  we  could  compute.  This  is  a  complete 
statistical  description  of  the  signal.  The  pre-processor  is  not  linear, 
however,  so  its  outputs  are  not  Gaussian.  Nevertheless,  the  joint  range- 
intensity  (or  Doppler-intensity)  probability  density  function  (pdf)  of  the 
output  random  variables  is  a  complete  statistical  description  of  those 
random  variablse.  This  chapter  derives  this  joint  density.  We  will 
perform  the  development  in  detail  for  a  ranging  radar;  results  for  the 
Doppler  and  range-Doppler  cases  are  handled  via  duality. 


A.  Signal  Model 

Figure  3.1  is  a  complex  envelope  model  of  the  ranging  radar  receiver 
through  the  pre-processor.  The  last  chapter  derived  an  equation  for 
j£_(  t  ).  For  a  ranging  radar,  the  radar  waveform  s^(  t  )  will  be  periodic 
with  period  T,  the  pixel  dwell  time: 


s.(  t  )  -  ^  J>(  t  -  nT  ), 


(  3.1  ) 


n  =  -  * 


where  J}(  t  )  is  some  elementary  waveform,  non-zero  only  in  the  interval 
[  -T/2,  +T/2  ].  The  filter  h_(  t  )  is  matched  to  the  transmitter  waveform 
b(  t  ): 


h(t)  =  b(T-t). 


(  3.2  ) 


The  next  block  is  a  square-law  envelope  detector.  In  many  real  systems 
this  block  would  be  a  linear  envelope  detector  because  the  linear  detector 
has  a  smaller  dynamic  range  at  its  output.  The  squared  magnitude  of  a 
complex  Gaussian  random  variable  gives  exponential  statistics,  whereas  the 
linear  magnitude  gives  Rayleigh  statistics  [16].  The  two  are  related  by  a 
simple  transformation,  but  the  exponential  statistics  are  easier  to  work 
with.  The  peak  occurs  at  the  same  point  in  either  case. 


Substituting  equation  3.1  for  js(  t  )  into  equation  2.3  for  j^(  t  )  gives: 

00 

JL(  t  )  ■  y  J~pI  I  dz  f  dp  _£_2(  p'  ,  Z,  t  -  z/c  ) 

V  0  i 

•  T(  p  ,  z,  t  -  z/c  )  _b(  t  -  nT  -  2z/c  )  +  _n(  t  )  .  (  3.3) 

The  matched  filter  output  is: 

a. (  t  )  =  x(  t  )  *  h.(  t  ) 

=  _x(  t  )  *  h.  (  t  )  +ji(t)*ji(t) 


=a(t)  +  a  .  (t), 

— x  —noise 


where 


a 

— x 


(  t  )  =  y  ->/  PT  /  dz  /  do  f 

4-*  V  0  ;  t  - 


t  -  T/2 


3T/2 


(  3.4  ) 


2  _ 

d  T  g  C  0,z,  T  -  z/c) 


T(  p  ,  z,  x  -  z/c  )  b(  t  -  nT  -  2z/c  )b(t  -  t  +  T  )  .  (  3.5  ) 


If  the  duration  of  b(  t  )  is  quite  short  compared  to  T,  T(  *  )  and  g  (  *  ) 
will  be  approximately  constant  over  the  period  of  time  b_(  t  )  is 
significant,  even  if  they  are  not  constant  over  the  entire  pixel  dwell 
time.  This  short-duration  pulse  condition  is  the  usual  route  to  obtaining 
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good  radar  range  resolution,  although  some  microwave  and  optical  radars 
use  long  duration  high  time-bandwidth  product  waveforms  which  undergo 
compression  later  in  the  receiver's  IF  system.  This  thesis  will  use  the 
short  pulse  assumption,  under  which  we  can  take  the  _T(  *  )  and  £  (  *  ) 
terms  out  of  the  time  integral,  so  that  a  (  t  )  is  approximately: 


a.x(  t  )  =  y  'sj  PT  I  dz  J  dp_£_2(  o  ,  z,  nT  +  z/c  ) 


•  T(  p ,  z,  nT  +  z/c  )  J 


t  -  T/2 
t  -  3T/2 


C  3.6  ) 


dx  b(  t  -  nT  -  2z/c  )  b  (  t  -  t  +  T  ). 


Meanwhile,  the  noise  component  is: 


a  .  (  t  )  -  n(  t  )  *  h(  t  ) 

-noise  — ^  — v 


rt  -  T/2  * 

=  /  d  t  n(  t  )  b  (  i 

J  t  -  3T/2 


-  t  +  T  )  . 


(  3.7  ) 

(  3.8  ) 


Since  ^_(  t  )  is  a  zero-mean  circulo-complex  Gaussian  random  process  and 
h(  t  )  is  a  linear  filter,  a(,  t  )  is  also  a  zero-mean  circulo-complex 
Gaussian  random  process.  Therefore,  its  magnitude  squared  process  will 
have  an  exponential  density  function. 

We  will  assume,  as  we  did  in  proposing  the  model  in  the  last  chapter,  that 
each  pixel  contains  only  one,  range-unresolved  reflector  element  (target 
or  background).  Again,  for  a  range  unresolved  target,  T(  ~  ,z,  t  )  = 

T(  p"  ,  L  ,  t  )  6  (  z  -  L>n  )  with  the  target  at  range  for  the  n-th 
pixel.  Now  the  filtered  signal  return  is: 


-  36  - 


a  (  t  ) 
— x 


=  T  y?~  J  dp  _g_2(  p  »  Ln,  nT  +  L Jc  )  T(  p  ,  Lfi,  nT  +  Ljc  ) 
n  *T 

,t  -  T/2  * 

*  /  dx  b(  x  -  nT  -  2L  /c  )  b  (  :  -  t  +  T  ) 

^  t  -  3T/2  n 


E  Sn  V  C  -  (n+1)T  '  2VC  > 


(  3.8  ) 


where  u  is  a  zero— mean  circulo-complex  Gaussian  random  variable  arising 

from  the  p -integral  and  R^C  x  )  is  the  autocorrelation  function  of 

_b(  t  ) ,  arising  from  the  x -integral.  Define  for  future  use,  t  n  =  (n+l)T 

+  2L  /c. 
n 


Finally,  we  divide  a_(  t  )  by  the  square  root  of  the  n(  t  )  white  noise 
spectral  density,  h  v  J  n  ,  to  get  a  normalized  signal: 

a'  (  t  )  =*  a.(  t  )  /-\J h  v  J  n 

=  ^/n  pt  /  hv  o  X  ^  Rbb(  t  -  Tn  )  +  ^  ) 

n 

*  a;  vc  c '  T- > c  y  (  3-‘ 


(  3.9  ) 


Here  v_(  t  )  is  the  normalized  noise  process  and  u^'  is  the  normalized 
target  return.  It  is  easy  to  show  these  have  moments: 


E[  u'  ]  -  0 


E[  u'  u* ,  ]  =  5  , 

— n  — n  nn 


E[  v(  t  )  ]  -  0 


E[  u'  u\  ]  -  0 
L  — n  — n 


E[  |v(  t  )|Z  ]  =  1 
E[  v(  t  )2  ]  =  0. 


(  3.10  ) 


v'.-> 
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Now  all  the  physics  of  the  transmitter  power,  LO  induced  shot  noise,  the 
beam  patterns,  the  target  reflectivity,  and  1/L  range  losses  are  in  the 
single,  measurable  parameter  CNR,  the  carrier-to-noise  ratio.  Physically, 
CNR  is  the  ratio  of  the  average  reflector  element  return  power  to  the 
average  LO  shot  noise  power  in  the  IF  bandwidth.  It  is  a  quantity  that  is 
well  known  in  microwave  and  optical  radar  theory.  For  the  monostatic 
radar,  resolved  speckle  target  case  of  interest,  CNR  obeys  the  radar 
equation  [17]: 


CNR  = 


h  ”  B 

o 


c  opt  ' 


(  3.11  ) 


where  B  is  the  bandwidth  of  the  IF  filter,  a  is  the  atmospheric 

extinction  coefficient,  --  q  is  the  radar  optical  efficiency,  and  c  is 

♦ 

the  target  reflectivity.  The  statement  E[  u  *u  ,  '  1=5  ,  also  implies 

all  pixels  are  independent.  This  is  true  because,  as  stated  earlier,  the 
scanning  assures  minimal  overlap  of  the  beam  between  two  pixels  and  the 
target  surface  has  a  very  short  correlation  distance.  Also,  since  the 
signal  and  detector  noise  are  from  different  physical  phenomena,  they  are 
independent:  E[u^'v(t)  ]  =  E[  a^' (  t  )  _v(  t  )  ]  =  0.  Because  the 
pixels  are  independent,  it  is  possible  to  examine  them  individually.  From 
now  on,  the  discussion  will  center  on  the  density  functions  for  a  single 
pixel,  so  we  will  arbitrarily  chosse  the  n  =  -1  pixel  and  drop  the  pixel 
subscript,  n: 


a'(  t  )  -  ^/CNR  u'  Rfeb(  t  -  T  )  +  v (  t  )  , 


T 

<  t  <  T  +  J  . 


(  3.12  ) 


This  equation  also  assumes  for  any  specific  time  t,  ^^bb  t  -  )  ~ 

(  t  -  T  )•  This  is  a  consequence  of  the  short  duration  of  the  radar 
waveform  Jj(  t  )  compared  to  the  dwell  time  T  and  the  assumption  of  a  range 
unresolved  target.  Put  together,  these  assumptions  assure  there  is  only 
one  return  per  dwell  time  and  this  return's  time  duration  is  much  less 
than  the  dwell  time.  Now  take  the  single  pixel  return  through  the 
squarer : 

w(  t  )  =  |a'(  t  )|  2 

=  CNR  |u'|2  |Rbb(  t  -  t  )|  2  +  |v(  t  )|  2  + 

2Re{  CNR  u'  Rbfa(  t  -  t  )  v*(  t  )  } .  (  3.13  ) 

The  intensity  process  w(  t  )  is  a  non-negative,  non-stationary  random 
process  whose  probability  density  function  for  any  single  time  t  is 
exponential. 

B.  Range-Intensity  Density  Function 

The  final  step  in  our  ranging  radar  pre-processor  is  to  peak  detect 
w(  t  ) ,  generating  the  two  random  variables  I  and  L,  the  intensity  and 
range,  the  *  meaning  an  estimated  or  measured  quantity.  We  desire  an 
expression  for  £(  x»  y  )»  the  joint  pdf  of  these  random  variables.  To 
derive  an  expression  for  this  pdf,  we  will  examine  w(  t  )  more  carefully 
and  construct  a  model  for  the  waveform. 

Suppose  we  wish  to  locate  a  target  return  time  within  the  dwell  time  of  T 
seconds.  Then  we  define  the  range  uncertainty  interval  as  =  cT/2.  If 
the  pixel  dwell  time  T  is  greater  than  the  time  corresponding  to  the  range 


uncertainty,  2L^/c,  we  will  of  course  attempt  to  locate  the  target  within 
the  smaller  time  interval.  At  the  risk  of  confusion,  we  will  use  the  same 
notation,  and  T,  in  either  case.  Figure  3.2  pictures  the  situation  for 
T  the  full  pixel  dwell  time  for  the  n  =  -1  case  we  introduced  in  equation 
3.12.  In  the  absence  of  noise,  _v(  t  )  =  0,  w(  t  )  has  a  peak  value  of  CNR' 
t_u *  |  at  t  =  T  =  t  ,  where  x  =  2L/c  is  the  time  delay  associated  with 
the  range  measurement  L,  as  sketched  in  figure  3.2(a)  for  an  arbitrary 
_b(  t  ).  If  _v(  t  )  £  0  but  the  CNR  is  very  large,  we  might  expect  w(  t  ) 
to  display  only  minor  corruption  from  the  noise,  as  in  figure  3.2(b),  so 
that  t  -  x  would  prevail.  However,  since  _u*  is  a  random  variable  and 
v(  t  )  is  a  random  process,  there  may  be  intervals  where,  even  at  high 
CNRs,  the  realization  of  ju_'  |  is  very  small  while  _v(  t  )  takes  on  a 
larger  value  somewhere  in  the  interval.  In  this  case,  the  peak  in  the 
interval  occurs  at  some  x  far  removed  from  x  ,  the  location  of  the  noise 
free  peak.  Figure  3.2(c)  depicts  this  situation.  The  situation  pictured 
in  3.2(c)  is  called  an  anomaly,  designated  by  the  event  {  x  ^  } .  The 
situation  in  figure  3.2(b)  is  called  no  anomaly,  designated  by  the  event 
(  x  q  } .  The  statistics  of  I  and  L  conditioned  on  ~x  q  and  x  ^  will  be 
very  different.  The  two  e'  ents  are  mutually  exclusive  and  collectively 
exhaustive,  so  Pr(  Hq  )  +  Pr(  )  =  1. 

The  existence  of  anomalies  is  not  unique  to  this  problem;  it  is  a  feature 
of  nonlinear  estimation  problems  in  general  [4],  Another,  related, 
feature  of  nonlinear  estimation  problems  is  the  existence  of  a  threshold 
effect.  In  classical  nonlinear  parameter  estimation,  the  estimator's 
performance  degrades  gracefully  as  the  CNR  falls  until  it  reaches  some 
critical  point,  the  threshold.  Beyond  the  threshold,  performance  degrades 


more  quickly  as  CNR  falls  [4J.  This  is  a  very  general  result  in  nonlinear 
estimation  theory.  For  single  pixel  measurements  of  speckle  targets, 
however,  the  threshold  is  not  apparent  [23].  For  this  case,  the  anomalies 
always  occur  and  we  are  effectively  always  below  threshold.  For  the  case 
of  multipixel  speckle  targets,  however,  there  is  a  threshold  effect  [23], 
so  we  have  taken  some  time  to  discuss  the  phenomenon  here. 


Both  the  signal  portion  of  w(  t  )  and  the  noise  portion,  _v(  t  ),  have 
widths  or  decorrelation  times  (see  equation  (3.7)  )  on  the  order  of  tr> 
the  width  of  the  pulse  _b(  t  ).  So,  w(  t  )  is  approximately  constant  over 
time  intervals  less  than  t  and  points  separated  by  greater  than  tr  are 
approximately  statistically  independent.  Consequently,  we  will  divide  the 
interval  T  into  an  integer  number  of  subintervals,  Q,  each  of  length 
T/Q  =  t  and  model  w(  t  )  as  a  sequence  of  Q  independent  exponentially 
distributed  random  variables  as  shown  in  figure  3.2(d)  [12],  The  random 
variables  all  have  mean  1  when  the  reflector  element  is  not  in  the  bin  and 
mean  CNR  +  1  when  the  reflector  element  is  in  the  bin.  Let's  reiterate 
the  Q  bins  need  not  cover  the  entire  pixel  dwell  time  of  T  seconds. 
Instead,  the  bins  may  cover  only  a  subinterval  of  T  in  which  the  user  is 
interested. 


At  this  point  we  can  overlay  the  range  and  subframe  models  from  the  last 
chapter  onto  the  bin  model.  From  figures  2.4  and  2.5  we  can  see  the 
background  bins  must  be  the  last  (farthest  range)  bins  and  the  target  must 
fall  in  one  of  the  remaining  (nearest  range)  bins.  We  will  divide  the  Q 
range  bins  into  potential  target  range  bins  and  Qg  background  range 
bins  as  in  figure  3.2(e).  For  a  target,  in  the  absence  of  noise, 
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1  <_  Q  <_  Q^.  For  the  background,  in  the  absence  of  noise,  Q^,+l  <.  Q  <.  Q» 
and  Q  =  QT  +  Qg. 

The  usual  approach  would  be  to  estimate  the  range  variable  x  by  first 
performing  a  global  maximization  by  choosing  the  bin  containing  the  peak 
and  then  performing  a  local  maximization  within  that  bin.  The  maximum 

A 

value  would  be  the  intensity  estimate  for  that  pixel,  I,  and  its  location 

/V  /S  A. 

will  be  the  range  estimate  for  that  pixel,  x  (or  L  =  cx  / 2).  In  general, 
it  is  this  global/local  approach  [4]  which  produces  the  threshold  effect 
mentioned  above.  As  the  CNR  drops,  the  global  estimator  will  make  errors 
(anomalies)  more  frequently  causing  the  local  estimator's  performance  to 
degrade  drastically.  However,  we  have  already  noted  that  for  the  speckle 
target  statistics,  the  anomalies  dominate  the  error  probabilities  at  all 
CNRs,  so  we  need  not  worry  about  the  threshold  effect  on  a  single  pixel. 

In  addition,  we  are  going  to  concern  ourselves  only  with  global  or  coarse 
range  estimates  —  estimates  of  only  the  bin  containing  the  target  with  no 
local  maximization  within  a  bin.  It  is  possible  to  do  the  fine  range 
maximization  and  derive  the  densities  for  the  measurements  [18],  but  they 
are  more  difficult  to  work  with.  Although  the  fine  range  results  indicate 
processor  performance  more  accurately,  they  provide  no  additional  insight 
into  the  processor.  Both  processors  exhibit  the  same  qualitative 
dependences  on  the  range  and  angular  resolution  parameters  [18].  From 
here  on,  we  will  refer  to  the  range  measurement  random  variable  by  its  bin 

A.  /N 

number  Q  =  2L/ctr* 

We  will  use  conditioning  to  compute  the  pdf  q(  x,  q  ): 


£*  *(  x,  q  ) 

I.Q 

A 

2*  *  (  x  1  q,  a  0  )  Pr(  Q  =  q  |  ct  Q  )  Pr(  a  0  )  + 

xlQ.a  o 

£»  .  (  x  |  q,  a  .  )  Pr(  Q  =  q  I  a  i  )  Pr(  a  !  )  •  (  3.14  ) 

I|Q.«  i 

The  joint  density  function  is  also  a  mixed  density  function.  The  random 

/v  a 

variable  I  is  continuous,  while  the  random  variable  Q  is  discrete. 
Consequently,  the  probability  function  for  Q  is  a  probability  mass 
function  rather  than  a  density  function. 


The  first  step  in  deriving  the  joint  pdf  is  to  compute  the  probability  of 
no  anomaly,  Pr(  aQ  )  =  1  -  Pr(  ).  Call  the  Q  exponentially 
distributed  random  variables  in  our  model  for  q  =  1  to  Q.  Assume,  with 

no  loss  in  generality,  the  target  is  in  bin  number  1.  Then  event 
occurs  if  Zq  Z^  for  all  q  £  1.  Use  the  bin  independence  and  directly 
integrate  over  the  exponential  densities  of  Z^  and  Z^  to  obtain 

Pr(  otQ  )  =  Pr(  Z2  <  Z1  ,...,  Zq  <_  Z1  ) 


-  EZJ  Pr(  Z2  <  Zj  |  Zl  )^~L  ] 


(  3.15  ) 


=  Ez  [  (  1  -  exp(  -  Zj  )  y  ] 


r(  Q  )  r( 


CNFT7 


mm 


r(  Q  + 


CMS  VI 


where  the  parameter  a  is  1  /  (  CNR  +  1  ).  This  parameter  recurs 

•a 

frequently  in  the  I  densities,  so  this  shorthand  is  convenient.  The 
reflector  element  CNR  can  be  for  the  target  or  background.  We  will 
designate  these  CNRs  as  CNRt  and  CNR^,  respectively,  when  we  wish  to  make 
a  distinction  and  write  their  corresponding  a  parameters  at  and  a^.  For 
large  CNR  [4]: 

Pr(  a0  )  =  1  -  a  (  log(  Q  )  -  1  /  2Q  +  0.577  )  .  (  3.16  ) 


Finding  the  density  (  x  |  q,  a  ^  )  requires  using  Bayes  law  and 

conditional  random  variables.  With  the  bin  model  we  have  introduced,  the 


conditioning  on  the  range  measurement  Q  corresponds  to  saying  the  peak 

A 

intensity  occurs  in  a  particular  bin,  say  the  q-th  bin,  Q  =  q. 

/v 

Conditioning  on  further  requires  Q  *  q  *  Qif  where  Q^,  i  =  t,  b,  is 

the  true  target  or  background  range  bin,  respectively.  Call  the  random 


variable  for  the  intensity  in  each  bin,  q,  Z^,  as  before.  The  desired 


density  is  the  joint  conditional  density  for  all  Q  intensity  values  after 

/\ 

integrating  out  all  but  the  one  where  Q  =  q: 


r<x>  rco  r  a 

£*, *  (  x  |  q,  a  i  )  “  J  dx,  ...J  dx  J 

l|Q»  a  l  1  0  1  J0  q  1  0 


dxq+l 


dxr 


0 


i  , »  (  x1 , . . . ,Xq  |  q,  a  j  )  • 


i  *  v  *ii 

z! ,  —  ,Zq|Q,  a  1 


(  3.17  ) 


But: 


\  7  .A  (  X1 . XQ  '  q’  a  1  )  * 

zi »  —  ,ZQ  |  Q.  at  x 


(  3.18  ) 


(  Pr(  ct  1|Z1«x1,...,Zq«Xq,(Hi)  Pr(Q=q|Zj=Xj, . . .  ,Zq=*Xq) 


vSb 

r 

SS 


8 


#  -E-Z  Z  'X1 

9  •  •  •  9  *■ 
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(X. y • • • f X.-.  )  ) 


Pr(  Q=q  |  ct  T  )  pr(  a  L  ) 


The  first  and  second  terms  in  the  numerator  are  both  indicator  functions: 
they  are  either  1  or  0  depending  on  the  values  of  their  conditioning 
variables.  Together  they  force  all  Q-l  integrals  to  have  upper  limits  of 
Xq.  Since  all  the  Z's  are  independent,  the  last  term  in  the  numerator  is: 


£z z0(  V*“*xq  >  -  %  (  xq  }  17  £zi  xq  } 

W  ^i  q  =  1  q 

q  *  Q-l 


(  3.19  ) 


where 


J>7  (  x  )  =  —  exp( - )  u(  x  )  ,  Z  .  =  < 

J  z.  z.  Ml 


CNR  +  1  for  j  =  Q. 


for  j  *  Q, 


The  first  term  in  the  denominator  is  just  1  /  (  Q  -  1  )  since  an  anomaly 
is  equally  likely  in  any  of  the  bins  Q  =#  Q^.  The  integration  is  easy 
because  of  the  indicator  functions,  yielding: 

£*  «  (  x  |  q,  a  .  )  =  5— ? - -  e  (  1  -  e  )  (  1  -  e  )x 

l|Q»  a  x  m  alj 


u(  X  )  (  1  -  5  )  . 


(  3.20  ) 


The  calculation  for  £j|q  a  (  x  |  q?  a  q  )  is  only  slightly  different: 


£•  »  (  x  |  q,  aQ  ) 

1 1 Q»  ex  0  U 


(  3.21  ) 


,r dx.  ...r dxal/  d xa+1...r 

JQ  1  J0  q_1  Jo  q+1  Jo 


. .  ./  dxQ  6 


Q  uqQi 


-  46  - 


Pr(  a0  I  Q*q,Z1=x1,...,ZQ=,XQ  )  2.  xx . xQ  |  q  ) 


\ . zq|q  'I' 


Pr(  a  0  ) 


Again,  Pr(  a  q  |  Q=q,  Zj=x^,...,  Zq=Xq  )  is  an  indicator  function  which 
forces  the  upper  limits  on  the  integrals  to  x^.  The  result  is: 

.  (  X  I  q,  a  0  )  -  Pr< 4  7  a  (  1  -  e'X  ><M  u(  *  >  5  ,Q  ' 

l|Q»ct  o  u  1 

(  3.22  ) 


Turning  to  the  probability  mass  functions  for  the  range,  under  event  a  . 
any  range  outside  the  correct  bin  is  equally  likely,  so: 


Pr(  Q  =  q  |  a  j  )  -  q-i-T  <  1  -  «  qQl  )  • 


(  3.23  ) 


where  is  again  the  noise  free  (true)  target  range  bin. 


The  probability  mass  function  Pr(  Q  =  q  |  a  q  )  is  by  far  the  most 
trivial.  It  is  simply  an  indicator  function: 


Pr(  Q  =  q  |  aQ  )  =  6  qQ_ 


(  3.24  ) 


Putting  all  the  pieces  together,  the  joint  density  is: 


£«  »(  x,  q  )  =  a  e  aX  (  1  -  e~x  )Q_1  u(  x  )  <5  0  + 

I,Q  4Vi 


(  3.25  ) 


5"X  (  1  -  e'aX  )  (  1  -  e"X  )Q"2  u(  x  )  (  1  -  5  „n  )  . 


vN/  .  vV  v  ,*.\* >>’«./ 


TO 
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The  marginal  densities  are 


,  \  -ax  /  i  -x  \Q~1  t  \ 

p»(  x  )  =  a  e  (1-e  )  u(x)  + 

I 


e_x  (  1  -  e-ax  )  (  1  -  e'x  )Q"2  u(  x  ) 

=  al  [  (  1  -  e-ax  )  (  1  -  e~x  )Q_1  ]  U(  x  )  , 

Pr(  “1  > 

Pr(  Q  =  q  )  =  Pr(  aQ  )  6  qQ_  +  —— —  (  1  -  5  qQ _  )  . 


(  3.26  ) 


(  3.27  ) 


C.  Range-Doppler  Duality 

The  previous  derivations  centered  on  intensity  and  range  measurements. 
This  section  will  demonstrate  a  duality  between  the  range  and  Doppler 
shift  measurements  that  will  make  it  unnecessary  to  always  repeat  range- 
intensity  analyses  for  Doppler-intensity  radars. 


The  previous  analysis,  up  to  equation  3.3  for  t  ) ,  does  not  change.  We 
will  pick  up  the  Doppler  analysis  from  that  point,  but  drop  the  n 
summation  index  and  let  n  =  -1  so  as  to  deal  with  only  one  pixel  (as  we 
did  in  the  range  analysis).  A  well  known  result  of  Fourier  transform 
theory  tells  us  that  in  order  to  obtain  good  frequency  resolution  we 
require  a  long  observation  time.  So,  for  a  Doppler  radar  the  duration  of 
K  t  )  will  generally  be  comparable  to  the  pixel  dwell  time.  Under  this 
condition,  almost  any  target  is  range  unresolved,  so  we  can  let: 


T(  o  ,  z,  t  -  z/c  )  =  T(  p  ,  t  -  z/c  )  6 (  z  -  (  L  +  v  t  )  )  (  3.28  ) 
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where  v is  the  target's  longitudinal  velocity  (only  the  longitudinal 
component  contributes  a  Doppler  shift).  For  simplicity,  assume  a 
transverse  velocity  of  zero.  Performing  the  z  integration  in  the  single 
pixel  reflector  element  return  integral  then  gives: 


L+v  t 


L+v  t 


x(  t  )  J  dp_§_(p,L  +  v2t,t - )  T(  p  ,  t - ) 


b(  t  - 


2(  L+v  t  ) 


)  +  n.  (  t  )  . 


(  3.29  ) 


If  we  perform  the  integrations  to  propagate  the  beam  to  the  target,  we 
find  under  the  far  field  approximation: 


JL  (  p  ,  z,  t  )  = 


_  !S  e2jkz  ejkl  P  f2/z 


2  rr  a(  z  ) 


exp[  - 


p  -  z  <f>  C  t  ) 
2  a  (  z  )2 


(  3.30  ) 


a(  z  ) 


»  d  . 


We  will  continue  to  use  the  assumptions  the  target  changes  little  over  a 
pixel  dwell  time,  and  that  there  is  negligible  radar  lag  angle.  We  shall 
also  assume  that  the  dwell  time  is  short  enough  that  only  a  relatively 
small  fractional  change  occurs  in  the  area  the  radar  illuminates.  Under 
these  conditions  we  have  that 


J_(p  ,z  +  vzt,t)  -  (  p  ,  z,  t  )  e 


jkv  t 


(  3.31  ) 


and 
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2  L  2jkv  t 


JL(  t  )  =  J?r  b(  t  (  1  -  ) - )  e 


/  dp  X.2(  p  .  z.  t  -  \  )  T(  7  ,  t  -  \  )  +  n(  t  ) 

Ap 


2  L  x  2JkV 


pT  u  b_(  t - —  )  e 


+  n  (  t  )  , 


(  3.32  ) 


where  we  used  2vz  <<  c.  Here  b_(  *  )  is  the  delayed  pulse,  u_  is  the 
speckle  induced  amplitude  variations,  and  the  exponential  is  the  Doppler 
shift. 


Now  let  us  impose  a  pre-processor  as  before  except  instead  of  matched 
filtering  jr.(  t  ),  perform  a  time-limited  (windowed)  Fourier  transform.  So 
the  output  of  this  transform  will  be: 


f  >  -  FV  hV  x(  t  )  ) 


u1  f  W(  t  )  b(  t  -  ~  ) 

*  CO 


2  L  ,  j2kvzt  “J27ft 


+  v(  f  )  , 


(  3.33  ) 


where  u/  is  \i  scaled  to  a  unit  variance  and  W(  t  )  is  the  window  function, 
W(  t  )  =  0  for  j 1 1  >  T/2.  The  noise  v(  f  )  is  a  unit  variance,  zero  mean, 
complex  Gaussian  random  process.  Continuing: 


a’(  f  )  =  -J CNR  u'  [  W(  f  -  f .  )  *  b(  f  -  f  )  ]  e 


-J2-t  (  f  -  fd  )t 


(  3.34  ) 


+  v  (  f  )  . 


* 

Note  that  if  we  choose  W(  t  )  so  that  W(  f  )  =  J)  (  -f  ) ,  the  convolution 
becomes  Rgg(  f  -  )  and  we  have  an  exact  analogy  with  the  ranging  radar 

result  in  equationn  3.18,  except  for  the  phase  term.  However,  we  can 
combine  the  phase  term  with  u/  without  changing  its  statistics.  So,  we 
have  demonstrated  an  exact  analogy  between  the  ranging  and  Doppler  radar 
signals  under  the  stated  conditions.  By  using  appropriate  scaling 
factors,  we  can  use  all  the  range  radar  results  in  Doppler  radar 

A  A  A 

analyses.  The  Doppler  Q  measurement  is  Q  =  F  /  f  where  the  random 
variable  F  is  the  measured  Doppler  shift  and  f  is  the  frequency 
resolution  or  Doppler  bin  width  in  Hz. 

Although  we  have  established  a  duality  between  the  IF  pre-processor  output 
signals  in  ranging  radars  and  Doppler  radar,  we  must  still  demonstrate  a 
duality  between  the  target  and  background  models  for  the  two  systems.  We 
will  still  assume  the  down-looking  radar  geometry  of  figure  2.3.  If  the 
radar  is  moving  parallel  to  the  ground  with  velocity  V,  we  can  use  simple 
trigonometry  to  determine  the  projection  of  v  onto  the  line  of  sight 
direction  to  the  ground  in  terms  of  the  radar's  pointing  angles,  T(  t  ). 
So,  we  can  assume  we  know  the  longitudinal  velocity  of  the  background 
relative  to  the  radar  for  each  pixel.  If  the  radar  is  stationary,  we  have 
the  special  case  V  =  0  and  all  the  background  pixel  velocities  are  0.  As 
long  as  the  target  moves  with  a  constant  longitudinal  velocity  with 
respect  to  the  ground,  its  velocity  will  fall  into  a  Doppler  bin  distinct 
from  the  background  bins.  However,  since  the  target  may  move  toward  or 
away  from  the  radar,  it  is  not  possible  to  place  the  Qg  background  Doppler 
bins  at  either  end  of  the  Q  total  bins  in  general.  Typically,  one  would 
arrange  the  frequency  uncertainty  interval,  F^,  the  dual  of  the  range 
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uncertainty  interval  Lu,  so  it  is  centered  around  the  background  Doppler 
bins  unless  we  have  a  priori  knowledge  of  the  sign  of  the  target's  Doppler 
shift. 

D.  Range-Range  and  Doppler  Dulaity 

In  the  previous  section,  we  developed  the  duality  between  range  radar  and 
Doppler  radar  by  putting  both  post-filtering  waveforms  into  the  form: 

a'(  t  )  -  -y/CNR"  u*  Rbb(  t  -  x  )  +  v (  t  )  (  3.35  ) 

with  appropriate  joint  statistics  for  _u'  and  _v(  t  ).  However,  for  the 
peak  detecting  pre-processor  we  will  deal  with,  this  is  a  bit  more 
stringent  than  we  really  require.  Although  certainly  the  statistics  must 
be  the  same,  the  function  Rbb(  t  )  need  not  be  an  autocorrelation 
function. 

For  a  receiver  measuring  both  range  and  Doppler,  the  maximum  likelihood 
processor  is  the  ambiguity  processor  [19]  which  first  extracts  the  time- 
frequency  correlation  function 

-x(  td’  fd  )  =  /  dt  x(  t  )  s  (  t  -  td  )  e  (  3.36  ) 

—  00 

prior  to  envelope  detection  and  peak  detection.  If  we  choose  a 

transmitter  waveform  s(  t  )  such  that  the  ambiguity  function 

2 

I  x  (  td*  fd  >1  •  where 

r+  00  * 

X (  V  fd  )  "  /  dt  s(  t  )  s  (  t  -  td  )  e 


—  00 


(  3.37  ) 


-A. 


•> 
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is  the  near  ideal  thumbtack,  as  Skolnik  [20]  calls  it,  shown  in  figure 


3.3,  then  the  pre-processor  input  is: 


a'(  t,  f  )  =  ^/CNR  u'  x  (  t,  f  )  +  v(  t,  f  ) 


(  3.38  ) 


where  v(  t,  f  )  is  a  zero  mean,  unit  variance,  complex  Gaussian  random 
process  arising  from  the  filtered  local  oscillator  shot  noise. 


Now  we  divide  the  t-f  plane  into  Q  range-Doppler  bins.  The  waveform 
s^(  t  )  sets  the  range  and  Doppler  resolutions,  through  x  (  t,  f  ),  and  so 
the  bin  sizes.  Simply  mapping  the  2-D  array  of  Q  bins  onto  the  1-D  array 
of  Q  range  bins  makes  the  range-Doppler  receiver  based  on  equation  3.38 
isomorphic  to  the  range  radar  receiver  based  on  the  coarse  range  version 
of  equation  3.12.  Notice  in  both  this  and  the  range  only  processor,  by 
using  only  coarse  range,  Doppler,  or  range  and  Doppler  data,  the  exact 
form  of  the  R^C  *  )  or  X  (  *  )  functions  are  not  very  important  as  long 
as  the  noise  free  peak  is  unambiguous.  If  we  had  used  fine  range  data  as 
well,  the  analogy  would  need  to  be  more  stringent  to  translate  the 
results.  The  exact  shape  of  the  correlation  functions  would  need  to  be 
the  same. 


This  establishes  the  duality  in  the  IF  signals,  but  we  need  to  establish 
the  duality  between  the  range-Doppler  model  and  the  range  model  as  well. 
This  is  very  easy  now  that  we  have  established  the  duality  in  the  Doppler 
radar  case.  Taking  a  composite  of  the  range  and  Doppler  radar  models  we 
have  already  developed,  we  can  segment  the  range-Doppler  bins  into  a  set 
of  Qg  large  range,  low  Doppler  shift  bins  where  the  background  returns 
must  fall  and  the  remaining  Q^,  bins  where  a  target  could  occur,  as  shown 
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Figure  3.5 

Thumbtack  Ambiguity  Function 


.  A 


in  figure  3.4. 


With  the  joint  and  marginal  densities  in  hand,  we  can  proceed  with 
processor  design  in  the  next  chapter.  Using  the  dualities  we  just 
derived,  we  only  need  do  this  design  once  for  the  range  receivers. 


IV.  BINARY  DETECTION  RECEIVERS 


A.  The  Neyman-Pearson  Criterion 

This  thesis  focuses  on  binary  detection  receivers  which  use  the  Neyman- 
Pearson  criterion  for  detection  [4].  The  Neyman-Pearson  approach  is 
useful  for  deciding  between  two  events  whose  a-priori  probabilities  are 
either  unknown  or,  as  in  this  case,  not  meaningful  because  they  are  not 
random  phenomena.  The  problem  here  is  to  choose  between  two  possible 
decisions:  target  present  (H^)  and  target  absent  (Hq)  based  on  an 
observed  data  vector  R.  Here  R  may  be  any  subset  of  the  measured 

A  /V 

quantities  {  I  ,  Q^:  1  £  m  _<  M,  1  <.  n  N  }  where  the  index  m  represents 
the  subframe  number  and  the  index  n  represents  the  pixel  number  within  a 
subframe.  The  Neyman-Pearson  criterion  constrains  the  probability  of 
false  alarm,  Pp  ■  Pr{  declare  |  Hq  } ,  to  be  less  than  or  equal  to  some 
specified  value  Pp*  and  then  minimizes  the  probability  of  miss,  PM  = 

Pr{  declare  Hq  |  H^  } ,  subject  to  the  constraint.  Performing  this 
constrained  minimization  results  in  a  likelihood  ratio  test: 


JES|hx(  ?  I  Hi  } 

^RIHq^  r  I  Ho  ) 


(  4.1  ) 


with  A  chosen  so  Pp  =  Pp’ . 


If  A(  r  )  is  not  a  continuous  random  variable,  Pp  will  not  be  a 
continuous  function  of  A  .  In  this  case,  if  the  desired  value  of  Pp  is 
not  one  of  the  discrete  values,  we  must  randomize  the  decision  rule  to 
obtain  the  desired  performance  [4],  The  reference  describes  the  procedure 
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for  a  simple  example.  We  will  need  to  use  this  procedure  in  chapter  6  to 
analyze  the  performance  of  the  range  processor  because  the  discrete  range 

A 

measurements,  {  Q  },  give  rise  to  a  discrete  likelihood  ratio  and  Pr. 

mn  r 

function. 


B.  Unknown  Parameters 

Equation  4.1  is  not  completely  general,  since  the  densities  often  depend 
not  only  on  the  data  R,  but  also  on  some  other  unknown,  non-random 
parameters.  If  it  is  impossible  to  express  the  likelihood  ratio  A(  r  ) 
independent  of  these  parameters,  we  must  first  determine  the  maximum 
likelihood  estimates  of  the  parameters  under  Hq  and  and  then  use  these 
estimated  values  in  the  densities  as  if  they  are  the  true  values.  This 
procedure  generates  the  generalized  likelihood  ratio  test  [4]: 


A  (  r  )  s 

°  max  { 


"S'11  ;  I  Hj,  A  )  }  *1 

- - - - — - - - 


fxl  *8|Hn,A(  r  I  V  »  )  ) 


(  4.2  ) 


where  A  is  the  unknown  parameter  vector. 


In  this  chapter  we  will  specifically  treat  intensity  and  ranging  radars 

and  later  treat  Doppler  radars  by  duality  arguments.  For  the  intensity 

and  range  radars  we  are  analyzing,  the  unknown  parameters  are  mQ,  the 

actual  subframe  which  contains  the  target,  ()  ,  the  actual  target  bin 

tran 

number  for  the  n-th  pixel  in  the  m-th  subframe,  CNRt  and  CNR^,  the  target 
and  background  carrier-to-noise  ratios,  respectively,  and  sometimes 
functions  of  these  parameters  (i.e.,  Pr(  ctgt),  the  probability  of  no 
target  anomaly,  a  function  of  CNRt). 


Introducing  the  generalized  likelihood  ratio  is  a  crucial  step  in  this 
development.  The  unknown  parameters  involve  the  multipixel  nature  of  the 
target.  Introducing  the  generalized  likelihood  ratio  in  combination  with 
the  models  of  chapter  2  is  the  procedure  which  makes  it  possible  to  derive 
optimal  or  quasi-optimal  processors  for  multipixel  targets.  Although  the 
generalized  likelihood  ratio  concept  is  not  new  [4],  applying  the  concept 
to  this  problem  is  new  and  this  is  what  allows  us  to  derive  these  results 
for  multipixel  processors. 


C.  Subframe  and  Frame  Statistics 

This  section  presents  the  general  method  we  will  use  to  derive  the  frame 
statistics  from  the  subframe  statistics  and  treat  the  target  true  subframe 
location,  mQ,  as  an  unknown  parameter. 


We  will  have  some  frame  level  measurement  vector,  which  is  a 

’ —  (s) 

composite  vector  of  M  subframe  measurement  vectors,  R^  ,  where  the 
superscripts  (f)  and  (s)  indicate  frame  and  subframe  level  measurements, 
respectively.  The  generalized  likelihood  ratio  test  with  unknown 
parameters  mQ  and  A  (for  any  remaining  unknown  parameters)  is: 


(f) 


max  {  I  H,  ,  A,  mA  )  )  H1 

A,  mQ  R^lHj.5,-  1 


,m_ 


> 


.max  {  £_(n  _ 

A,  mQ  RU;|H0,A,mo 


(  ?(f)|  H0,  A,  mQ  )  }  H 


X  ,  (  4.3  ) 


0 


Since  the  subframes  do  not  overlap,  the  subframe  measurements  are 
statistically  independent  and  those  subframes  not  containing  the  target 
have  the  same  statistics  as  if  there  is  no  target  in  the  frame: 


-  59  - 


“J  ^S)|H1,5 


(  r(ms)  |  H  A  ) 
m_  L 


n  w.  *  ?“s)  1  h°'  5  j 

m  =  1  m  |H0,A 


m  *  m 


i“  H  V)|Vi(  ?»s)|ho-5> 


.  (  4.4  ) 


Since  the  denominator  in  equation  4.4  is  independent  of  mo>  we  can  drop 

the  m  maximization  in  the  denominator  and  take  the  denominator  inside  the 
o 

m^  maximization  in  the  numerator  (being  careful  to  keep  the  A 


maximizations  separate) : 


T  ^  '  n*’  1  ’ 


n 

•  n  V)|H  A(  ?"S>  1  "O’  5  > 
m  =  1  m  |H0,A 


m  ^  m 


max  <  1  ["  p  /  (  r^l  Hn,  A  ) 

A  li,  R(s)|h„,5  “  0 

m  =  1  m  1  0 


>.(  4.5  ) 


Often  the  unknown  parameter  vector,  A,  contains  elements  which  affect  the 
densities  under  or  Hq,  but  not  both.  Examples  are  CNRt  which  affects 
the  density  under  but  not  Hq,  and  CNR^  which  affects  the  density  under 
Hq  but  not  H^.  Let's  call  such  parameters  "hypothesis  separable" 
parameters.  Then  the  vector  A  becomes  (  Aq,  A^  )  where  the  subscripts 


ksd 
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indicate  which  hypothesis  the  parameters  affect.  If  the  parameters  are 
hypothesis  separable,  we  can  rewrite  the  part  of  equation  4.5  inside  the 


m  maximization  as: 
o 


max 

A, 


M 


max 

A, 


0  I 


n  £R(s)iHn,50( 

m  =  1  m  1  O’  0 


m  1  O’  0  J 


m  Am 


M 


••  (  4.6  ) 


ti  n  i ( ?”S)|  H°’  5° >! 

A0  m  =  1  m  ,H0’A0 


The  difference  between  the  Aq  maximizations  in  the  numerator  and 


denominator  is  one  subframe  term.  If  M  is  fairly  large,  both  numerator 
and  denominator  will  give  nearly  the  same  estimates  for  the  parameter 


vector  Aq.  If  we  make  the  estimate  in  the  denominator  and  use  the 


resulting  vector  An  in  the  numerator,  we  will  introduce  only  a  small 


error,  but  all  M  -  1  terms  where  m  £  m  will  cancel  in  the  numerator  and 

o 


denominator  leaving  for  equation  4.6 
max  { 


nax  i  p  ,  * 


T  (  ?iS)  I  H  Aj  )  ) 

0 


l»0-5 

o 


( ;i9)  >  ,  ( 4.7  > 


p  .  .  A  (  rK3f  |  Hn,  An  ) 

d(3)  lu  7  0  0 

0 


gn 


where  Aq  is  the  estimate  of  vector  Aq  derived  from  measuremnets  over  the 


entire  frame.  Notice  we  have  introduced  a  shorthand  notation  for  the 


(s) 


subframe  generalized  likelihood  ratio,  A  '  .  Now  we  may  write: 


-  61  - 


A(f)  *  max  {  A(s) 
8  m  8* 


(  4.8  ) 


Now  it  is  clear  the  generalized  likelihood  ratio  for  the  frame  is  simply 
the  maximum  of  the  generalized  likelihood  ratios  for  the  M  subframes. 

This  fact  makes  it  quite  easy  to  find  the  frame  level  probabilities  of 

false  alarm  and  miss,  Pp^  and  from  the  subframe  level 

( s )  (s) 

probabilities  of  false  alarm  and  miss,  Pp  and  .  The  procedure  is 


as  follows: 


=  Pr{  /^f)(  r(f)  )  >  A  |  HQ  } 


(  4.9  ) 


=  1  -  Pr{  <  A  |  Hq  } 


=  1  -  Pr{  max{  Al„'  }  <  X  |  HQ  } 
m  ® 


We  will  treat  all  M  of  the  Agm  '  as  approximately  independent  and,  under 

Hq,  identically  distributed.  They  are  not  exactly  inndependent  because, 

—  ( s )  — * 

although  the  measurements  Rm  are  independent,  the  Aq  value  used  in  the 
denominator  is  derived  from  a  maximization  over  measurements  for  the 
entire  frame.  If  the  maximum  is  less  than  some  specific  value,  all  M  must 
be  less  than  that  value,  so: 


-  Pr(  <  *  1  Ho  >M 

-  [  1  -  Pr{  a£5  >  A  I  H0  )  ]M 

-  [  1  -  P<9>  ]M  . 


(  4.10  ) 


Here  Pp^s'  has  an  intuitively  meaningful  definition.  It  is  the 
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probability  of  false  alarm  if  we  consider  just  a  single  subframe.  A 


nearly  identical  analysis  for  P. 


gives: 


PMf>  *  PHS>  [  1  -  PFS>  ^ 


again  assuming  independent  A 


/  M-l  v 
(  ncr  ) 


-  p^  r  i  i 

"  M  L  1  “  J 


(  4.11  ) 


where 


PM3)  -  Pr‘  <  x  I  Hx  )  • 


The  physical  interpretation  of  these  equations  is  as  follows.  A  miss 

(o') 

-  j c - ~c  «.u_  —  i-f- - '  '  clears  the  threshold 


occurs  if  none  of  the  subframe  statistics,  A„ 

gra 

X  .  A  detection  occurs  if  at  least  one  of  the  A 


statistics  clears 


the  threshold,  even  if  the  subframe  whose  statistic  clears  the  threshold 

does  not  actually  contain  the  target.  But  the  probability  of  a  detection 

on  the  wrong  subframe  turns  out  to  be  quite  small;  this  event  is 

(s) 

essentially  a  false  alarm  and  its  probability  is  bounded  by  Pp  .  Notice 

(s) 

this  argument  depends  only  on  Pp  being  relatively  small.  It  does  not 
depend  on  the  CNRs  or  any  other  parameters.  Further,  from  equation  4.10, 
if  M  is  large,  Pp^  may  be  small  even  if  Pp^  is  not. 


i  g  i 

Notice  since  adequate  performance  generally  requires  Pp  to  be  fairly 
small,  we  can  usually  make  the  approximations: 

p£f)  *  M  p£s) 


^  ’j-  i  -v-  -  v  >-  ---V- 


and 


(  4.12  ) 


p(s) 
*M  * 


In  words,  equation  4.12  says  for  a  fixed  target  size,  N,  Pp^  rises 
linearly  with  an  increase  in  the  angular  search  area  and  is 

approximately  independent  of  the  angular  search  area. 


D.  Multiple  Targets 

Suppose  we  wish  to  consider  the  possibility  of  multiple  targets  in  the 
frame.  We  will  continue  to  assume  targets  align  with  subframe 
boundaries.  We  will  also  assume  the  number  of  subframes  containing 
targets,  M' ,  is  a  small  fraction  of  M.  There  are  several  possible  ways  to 
analyze  this  problem  and  define  the  corresponding  false  alarm  and 
detection  (or  miss)  events,  each  of  which  may  make  sense  for  a  particular 
problem.  We  will  form  the  generalized  likelihood  ratio  for  each  subframe, 
perform  a  threshold  test  on  each  subframe,  and  declare  target  present  or 
target  absent  for  each  subframe  independent  of  all  other  subframes.  That 
is,  perform  no  maximization  over  the  subframe  generalized  likelihood 
ratios  before  the  threshold  test.  A  false  alarm  occurs  any  time  we 
declare  a  target  in  an  empty  subframe,  even  if  there  is  a  target  elsewhere 
in  the  frame.  A  miss  occurs  anytime  we  declare  target  absent  in  a 
subframe  containing  a  target  even  if  we  have  correctly  detected  targets  in 
other  subframes. 


With  this  processor,  the  probability  of  miss  for  any  target  is  just  the 
( s) 

value  already  computed.  Similarly,  if  we  compute  a  frame  level 

probability  of  false  alarm  Py^^'  it 


xs 


- 1  -  f  1  - 


Pr{  A (S0l  >  A 

gm  — 


Hq  j  r 


(  4.13  ) 


since  M' ,  the  number  of  subframes  containing  targets,  «  M  by  assumption. 


The  point  of  all  this  discussion  is  whatever  results  we  derive  for  the 
single  target  in  a  frame  are  also  useful  for  discussing  performance  with 
multiple  targets  in  the  frame.  We  need  make  only  minor  corrections  to  the 
single  target  performance  numbers  to  make  quantitative  statements  about 
multiple  targets  in  a  frame. 


E.  Measurement  and  Target  Models 

We  will  retain  the  pre-processor  introduced  in  chapter  2  and  continue 
working  with  only  coarse  range  or  Doppler  data  as  in  chapter  3.  In  this 
case,  we  already  have  the  density  functions  we  need  from  chapter  3.  Now 
we  will  use  these  density  functions  in  the  generalized  likelihood  ratio 
tests  from  above  to  derive  the  processors  for  receivers  measuring  range, 
Doppler,  intensity,  or  some  combination  of  these  signal  dimensions.  Since 
we  have  already  established  a  range-Doppler  duality  and  a  range  and 
Doppler  analogy  with  range  only  measurements,  the  joint  range-intensity 
measurement  receiver  is  our  most  general  case.  We  will  do  this  case  and 
then  the  range  only  and  intensity  only  cases. 


All  these  derivations  use  the  target  model  introduced  in  chapter  2.  In 
that  model  the  target  is  range  unresolved  at  an  unknown  range.  It  is 
vertical  against  a  background  sloping  away  at  a  known  range.  Neither 
target  nor  background  CNRs  are  known.  Targets  align  with  subframe 
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boundaries  and  we  wish  to  determine  whether  or  not  a  target  exists  in  a 
frame. 


We  will  label  each  pixel  with  a  subscript  for  subframe  index,  m. 


1  m  _<  M,  and  a  pixel  index,  n,  1  <.  n  jC  N,  as  we  did  at  the  start  of  this 


chapter.  For  each  pixel  in  a  subframe,  we  measure  the  intensity,  1^  , 

at  the  known  background  range,  Qg^.  This  gives  rise  to  an  exponential 

A 

density  function  for  1^  .  For  each  pixel  we  also  measure  the  intensity 

in  all  non-background  range  bins  and  choose  the  maximum  of  those 


intensities.  Call  this  maximum  I.  and  call  the  bin  it  occurs  in  (3  , 

tmn  tmn 

the  presumed  target  range.  This  gives  rise  to  the  statistics  we  derived 


in  the  last  chapter  for  I  and  (T  .  Of  course,  I,  is  independent  of 
r  tmn  tmn  bmn  v 

A  A  A  A 

both  I  and  ()  ,  though  I.  and  Q.  are  dependent  random  variables  in 

tmn  tmn  6  tmn  xtmn  v 

general.  This  is  the  most  general  form  of  the  model. 


F.  Joint  Range-Intensity  Processor 

If  there  is  no  target  present  for  a  particular  pixel,  the  system  CNR  will 
be  CNRg  and  the  effective  CNRt  is  zero.  There  are  Q  range  bins  and,  using 
the  notation  introduced  in  chapter  3,  one  of  Qg  of  them  is  the  background 
range  bin,  so  we  do  the  maximization  over  only  Q^,  «*  Q  -  Qg  bins  to  find 


■^tmn  anc*  ^tran  ^or  eac^  pixel*  We  measure  lgmn  at  the  known  background 
range  bin.  The  densities  for  the  resulting  measurements  are: 


I.  Hn 
bmn  0 


<  I  Ho  ) 


abXbmn  , 

ab e  uC  xbm  >  • 


(  x 

I.  ,Q.  |Hft 
tmn  tmn1  0 


tmn’ 


xmn 


Hq  >  - 


(  4.14  ) 


m 
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-x.  Qm-1  — x. 

,  .  tom  N  1  tmn 

(  1  -  e  )  e 


u(  *,mn  )  ,  1  <  qm  <  Qt 


The  background  and  target  random  variables  are  conditionally  independent, 


so  we  simply  multiply  the  two  densities  to  find  the  joint  density  for  the 


measurement  vector  (  1^,  ^tmn'  ^tmn  ^  on  eac^  Pixe^  under  hypothesis  Hr 


If  there  is  a  target  for  a  particular  pixel,  the  system  CNR  is  CNRt  and 


the  effective  CNR^  is  zero.  The  resulting  densities  are: 


1  “b™  I  H1  >  ■  e  u<  xbmn  >  • 

bmn'  1 


p»  »  (  xk  ,  q  H.  )  = 

■  v  tmn  ^mn  1 


Itmn,QtranlHl 


(  4.15  ) 


atXtmn  ,  <  Xtmn  1  ,  ,  ~ 

at  e  (  1  '  e  >  u(  xtmn  >  6 


q  Q  + 
4mn4t 


-a.x.  -x.  Qrp-2  — x. 

,  .  t  tmn  .  ,  .  tmn  .^T  tmn  ,  \  t  ^  &  \ 

(  1  -  e  )  (  1  -  e  )  e  u(  x  )  (  1  -  6  )  , 

Hmnxt 


*  —  ^mn  —  • 


The  background  and  target  random  variables  are  still  conditionally 


independent  under  hypothesis  H^, 


Since  all  pixels  are  independent,  the  frame  density  is  just  the  product  of 


the  M  subframe  densities  and  each  subframe  density  is  the  product  of  N 


pixel  densities  from  equation  4.14  and  4.15.  Applying  the  general  results 


of  part  C  above,  the  frame  generalized  likelihood  ratio  is 


A^  =»  max{  A  ^  }  , 

g  gm 


rv 


. 

mu 

i 

L  L. 

<  ‘bin  1  H1  >  ! 

bmn|nl 

(  4.16  ) 

n  =  1  tin  tan1  0 


Pa  »  (  x,  I  H„,  a,  ) 

t  |H  „  bmn  0  b 

bmn^  O’  b 


The  estimate  a^  is  from  the  maximization  over  an  entire  frame,  i.e., 
M  N 


=<  n  n  - 


ab  m  =  1  n  =  1  ^bnur^0‘ab 


(  Xbmn  I  H0’  ab  )  * 


occurs  at  a^  =  a^.  This  gives 


M  N 


A  Z  l  h 

m  =  1  n  =  1 


(  4.17  ) 


The  true  value  of  the  parameter  a^  is  the  true  mean  of  the  sampled  random 
variable,  1^  .  The  estimate  of  the  parameter  a^  is  exactly  what  we  would 

expect  from  intuitive  notions  and  the  law  of  large  numbers:  it  is  the 
sample  mean  for  measurements  of  1^  over  the  entire  frame. 


Now  we  substitute  the  densities  in  equations  4.14  and  4.15  into  equation 
4.16,  cancel  terms  in  common  in  the  numerator  and  denominator,  take  the 

f  a  \ 

logarithm  of  Agfflv  to  change  the  products  to  summations  and  use  the 
monotonicity  of  the  logarithm  to  justify  bringing  the  logarithms  inside 
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the  maximizations.  The  result  is: 


N 


log  A ^ s ^  =  max(  Y"*  log  A^(  x.  ,  q  ,  x.  )  }  ,  (  4.18  ) 

6  gm  /  .  6  v  turn*  Hmn’  bmn  J  1  ’  v 

V^t  n  =*  1 


where 


log  A^(  x  ,  q  ,  x.  )  =  A  ^  = 

6  tmn  Hmn  bmn  mn 


log 


rP*  *  *  (  xt  ,  q  ,  x,  I  H, ,  a.  ,  Q.  )  -i 

1  -z-j  n  T  j  tt  n  tmn  mn  bmn  1  1  t  t  | 

tmn’^tran*  bmn'“l ,at’^t 


I.  ,Q.  ,1,  Hn,a, 

tmn  tmn’  bmn  0  b 


(  xt  ,  q  ,  x,  I  Hn,  a,  ) 
tmn  mn  bmn  '  0  b 


[  logC  a  )  +  (  1  -  a  )  x  ]  S  q  + 

»  t 


(  4.19  ) 


[  log(  1  -  e  atXtmn  )  -  log(  1  -  e  Xtmn  )  ]  C  1  -  5a  0  }  " 

^mn't 


[  log(  ab  )  +  (  1  -  ab  )  xbffin  ]  . 


The  (p)  superscript  indicates  a  pixel  level  likelihood  ratio. 


Consider  a  very  large  CNRt  (say,  >60  dB  )  and  Itmn  CNRt,  which  then 
occurs  with  high  probability.  The  logarithm  of  A  is  then 

approximately: 

A^p)  =  -  [  log(  a,  )  +  (  1  -  a,  )  x,  ]  + 

mn  L  °  b  b  bmn  J 

(  4.20  ) 

[  log(  a  )  +  (  1  -  a  )  x  ]  5  q  • 

^ran  t 

Since  we  are  interested  in  receivers  we  can  analyze  rather  than  strictly 
optimal  receivers,  we  will  use  this  approximation  for  all  CNRts,  even 

-•v 

small  ones  (CNRt  >10  dB)  where  the  large  CNRt  and  Itmn  assumptions  are  not 
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valid. 


Now  we  perform  the  maximizations  over  Qt  and  at  in  equation  4.19  so  the 
subframe  log  likelihood  ratio  is: 


11 

lQg  ^  (gm  ma^(  £  108 


V^t  n  =  1 


A<P>  } 

mn 


- N  [  r  (  lo»(  *tm  >  +  —  -  1  >  - 

atm 


log(  ab  )  +  (  1  -  ab  )  Tbm  ]  . 


(  4.21  ) 


where 


y  5  . 

i  q  Q. 

n  =  1  Mmn  tm 


’  bm  ~  N 

n  =  1 


x.  6 

turn 


n  =*  1 


Q  Q 

mntm 


)  «  * 

^ — '  q  Q 

n  =  1  mn  tm 


M  N 


■  ab s  o  L  Z 


m  =  1  n  =  1 


bmn 


Physically,  j  is  the  number  of  times  the  peak  detector  selected  bin  Qtm 

as  the  maximum  bin  in  subframe  m.  We  will  coin  the  term  "hits":  i  is 

Jm 

the  number  of  hits  in  bin  Q.  .  Clearly,  j  is  an  integer  in  the  interval 

tm  Jm 

[  1,  N  ]  with  larger  values  more  probable  for  high  CNRt  values. 

A 

Physically,  afc  is  the  inverse  of  the  average  intensity  for  the  hits  in 
the  presumed  target  bin  Q  .  It  is  the  same  form  as  the  estimate  ab  (for 
a  single  subframe)  except  the  uncertainty  in  target  location  has  reduced 
the  number  of  intensities  to  average  over  from  N  to  j  .  The  quantity  Ibnj 


v  “^vyvyv/. 


s  v 


is  just  the  average  background  intensity  for  subframe  m. 

A 

The  only  thing  left  undefined  here  is  the  procedure  to  determine  Q  . 
Mathematically,  the  procedure  comes  from  the  at  and  Q  maximization  terms 

/V 

in  equation  4.21.  Physically,  the  mathematics  tells  us  bin  Q  is  the  bin 
which  has  the  largest  average  intensity  for  measurements  over  the  entire 
frame.  Figure  4.1  demonstrates  the  procedure  for  Q^,  =  6  and  N  =  5.  For 

/V  A 

each  pixel  we  record  Q.  and  I  .  For  each  bin  we  sum  the  intensities 
v  ^tmn  tmn 

of  all  the  hits  in  that  bin  number  over  the  entire  subframe  and  divide  by 
the  number  of  hits  in  that  bin  number  to  find  the  average  intensity  for 
that  bin.  We  compute  this  average  for  all  Qrp  potential  target  range  bins 

A 

and  select  the  maximum.  The  bin  where  this  maximum  occurs  is  Q  ,  the 

tm 

number  of  hits  in  that  bin  is  j  ,  and  the  inverse  of  the  average  for  that 

A 

bin  is  a  .  If  there  is  a  target  in  the  subframe  and  CNR  is  significant 

LIU  L 

(>  10  dB),  there  is  a  high  probability  Qtm  will  be  the  bin  with  the  most 
number  of  hits,  too.  If  CNRt  is  very  small  or  there  is  no  target  in  the 
subframe,  there  may  be,  for  instance,  one  bin  with  two  hits  and  a  smaller 
average  intensity  than  some  other  bin  with  only  one  hit. 

Equation  4.21  has  a  physically  intuitive  form.  It  is  essentially  the 
difference  between  the  average  intensity  at  the  presumed  target  range  and 
the  average  intensity  at  the  known  background  range.  Under  H^,  j 
should  be  reasonably  large  (  >N/2  )  and  the  a  dependent  terms  are 
approximately  CNRt  -  log(  CNRt  )  which  is  >  1  with  high  probability. 
Meanwhile,  the  background  terms  will  be  close  to  zero  with  high 
probability.  This  gives  a  large  log  likelihood  ratio.  For  hypothesis  Hq 
the  target  terms  will  be  small  and  the  background  terms  large,  so  the  log 


6  Pixel  #: 


I  =  1.13 
tmn 


X  =  "hit' 


ver  N  =  5  pixels/subframe 


a .  =  0.397 

tin 


likelihood  ratio  will  be  small.  If  MN  is  large,  we  can  treat  a^  as  an 
essentially  exact  estimate  and  thus  make  1^  the  only  background  dependent 
term.  This  simplifies  the  expression  some  since  we  can  add,  subtract,  and 
multiply  positive  constants  to  a  Neyman-Pearson  test  ratio  without 
changing  the  performance  (the  as  yet  undetermined  threshold  shifts, 
however).  An  equivalent  processor  is: 

~  /*■  <  ^  > +  —  ~ 1  >  -  Tb»  •  <4-22> 

1  -  ab  }  atm 

Here  I.  is  a  chi-squared  random  variable  with  2N  degrees  of  freedom  since 
it  is  the  sum  of  N  exponential  random  variables  [16].  However,  a is  a 
very  complicated  random  variable  and  enters  the  expression  in  a  very  non¬ 
linear  fashion.  Random  variable  j  is  only  slightly  less  difficult  to 
generate  a  probability  mass  function  for.  This  makes  it  very  difficult  to 
handle  the  performance  even  in  approximate  form.  We  could  always 
calculate  the  processors  performance  by  simulation  methods,  but  we  won't 
take  this  approach  either.  Instead,  we  will  examine  two  simpler 
processors  which  measure  range  or  intensity  only.  We  will  find  their 
performances  are  easier  to  handle  analytically  or  in  approximate  form.  We 
should  keep  in  mind,  however,  this  joint  range  and  intensity  processor 
will  out-perform  the  intensity  only  and  the  range  only  processors  or  any 
ad  hoc  combination  of  the  two  processors. 

G.  Intensity  Only  Processor 

In  this  section  we  will  investigate  a  receiver  which  measures  only  the 
intensity  of  the  pre-processor  output,  I  .  We  peak  detect  over  all  Q 
bins,  but  throw  away  and  we  don't  use  the  information  about  the  range 


.*» .%  .v 
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£l 


to  background.  We  use  only  the  measured  peak  intensity  I  .  We  will 
begin  the  intensity  only  processor  analysis  by  examining  the  processor  for 
a  single  pixel  measurement: 


max{  2*  (  x  j  H. ,  a  )  } 

T  I  TT  _  1  L 


A(P}(  x  ) 


at  I'Hl»at 


max{  (  x  \  H  a  )  ) 

ab  I'H0,ab 

,£*,  (  X  |  H.,  a  ) 

max  - 

at’ab  (  x  |  H0.  ab  ) 

1 1 M0  ’  ab 


(  4.23  ) 


A 

If  we  make  the  large  1^,  large  CNRt  and  CNRb  assumptions,  as  we  did  for 
the  joint  processor  earlier,  the  conditional  density  for  either  hypothesis 
reduces  to: 


(  x  |  H  ,  a  ) 
I  ,a .  1  J 


(  1  -  e  *  )  e  3jX  u(  x  ) 


(  4.24  ) 


for  i  =  0,  1,  and  j  =  b  for  i  =  0,  and  j  =  t  for  i  =  1 .  Then  the  log 
likelihood  ratio  is: 


log  A  £  (  x  )  =  max  (  (  ab  -  a  )  x  }  . 

Vab 


(  4.25  ) 


Since  the  ab  and  at  values  form  a  single  multiplicative  constant,  they  are 
irrelevant  to  a  Neyman-Pearson  processor  _if_  we  know  the  sign  of  ab  -  at« 
Then  the  generalized  likelihood  ratio  is  also  a  uniformly  most  powerful 


(UMP)  processor  [4]  and  is  simply: 


m 


log  X  )  =  X  . 

O 


(  4.26  ) 


Using  the  general  results  from  earlier  we  now  find  the  frame  level 
processor  for  N  pixel  subframes  is  simply: 


"f1  E  108  A 'gin  1 


n  =  1 


=  max{ 


y  x  } 

/  .  mn 


(  4.27  ) 


n  =  1 


max{  log  )  . 

m  ® 


where  it  is  still  critical  that  we  know  the  sign  of  ab  -  afc.  In  general, 

however,  we  don't  know  the  sign  of  a^  -  at<  Indeed,  we  may  have  either 

CNRb  >  CNRt  (called  negative  contrast),  or  CNRfc  >  CNRfe  (called  positive 

contrast),  or  CNRb  =  CNRt  (no  contrast  or  zero  contrast).  Nevertheless, 

in  keeping  with  the  idea  of  working  with  processors  we  can  analyze,  we 

C S )  (s)  r“ A 

will  continue  to  use  log  A  _ =  R  =  )  I _  as  our  subframe 

°  gm  m  mn 

statistic,  and  derive  a  way  to  handle  the  unknown  contrast  ratio. 


Since  we  will  generally  be  concerned  with  N  10  or  more  pixels  per 
subframe,  we  have  a  sum  of  a  "large"  number  of  independent  random 

(s) 

variables.  This  gives  us  a  conditional  Gaussian  density  for  R^  ,  via 

(s') 

the  Central  Limit  Theorem.  The  conditional  moments  for  R  are: 

m 

E[  R^S)  |  HjS)  ]  -  NMj,  var[  R^s)  |  H^s)  ]■»  oj  (  4.28  ) 


for  j  =  t,  b,  where  y  .  and  a  .  are  the  conditional  pixel  mean  and 

J  (s)  (s) 

variance,  respectively,  and  H  denotes  the  target  subframe,  and 

denotes  a  background  subframe.  After  some  straightforward  but  tedious 

calculations  we  get  the  conditional  moment  values: 
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(  4.29  ) 


where  ^ (  x  )  and  y 1 (  x  )  are  the  first  and  second  logarithmic 
derivatives  of  the  gamma  function,  p(  x  ).  The  equations  are  valid  for 
target  and  background  with  appropriate  subscripts  on  the  variables  u ,  a 
a,  CNR,  and  ctg. 


Since  the  subframes  are  independent,  the  frame  processor  falls  into  the 
same  general  form  we  derived  in  parts  B  and  C  above: 


"a,!(  K (3) 


“t  R; 


(  ri3’  I  Hj,  Ut  )  K 
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=  max 
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M3>IVb'  “ 


<  1  «n.  *  k  ) 


(  4.30  ) 


We  will  do  the  maximization  in  the  numerator  by  differentiating  the 
density  with  respect  to  u  directly  and  setting  the  derivative  to  zero: 


The  density  for  R^8'  ■  is  conditionally  Gaussian  with  mean  N 

2  2  n  2 

and  variance  N  afc  •  However,  a  is  a  function  of  yt  via  the  equations 
in  4.29.  This  dependence  between  p  and  at  complicates  the  formula  for 

A 

u  .  With  this  functional  dependence,  the  estimate  is  given  implicitly 


by  the  equation: 


utm  tm 


-  r  1  /  1 

atm  L  2  o’  +Vl  +  4a 


(  4.32  ) 


where  the  value  of  a  is  the  value  of  the  equation  for  at  evaluated  at 

A  -A 

Vi t  =  y tm  and  a'tm  is  the  derivative  of  at  with  respect  to  yfc 

/A. 

evaluated  at  y.  •  If  a*  was  not  a  function  of  y.  (or  CNR.),  the 
tm  t  t  t 

variance  dependent  terms  in  equation  4.32  would  not  be  present  and  the 

estimate  would  be  y_  =  y„  3  )  I  ,  the  sample  mean.  At  very  high 

tm  tm  *— 1  mn 

n 

CNRt  values,  at  -  ut»  so  the  estimate  reduces  to: 


=  0.618  y 


(  4.33  ) 


For  small  CNRt  values,  at'  and  at  approach  constant  values  and  u 
approaches  the  sample  mean  minus  a  constant  bias.  In  keeping  with  our  aim 
to  use  simple  quasi-optimal  processors,  we  will  use  the  approximation 

A 

y  *  i^tm  3Ince  is  a  reasonable  approximation  and  has  the  advantage 
of  being  independent  of  CNRt> 


To  obtain  y^,  we  perform  a  similar  differentiation  except  we  use  the 
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statistics  for  the  entire  frame:  , . . .  ,Rj^S^  )*  where  the 

(s) 

Rm  random  variables  are  conditionally  independent  jointly  Gaussian 

2  2 

random  variables  with  identical  means,  N  and  variances,  N  . 

Setting  the  derivative  with  respect  to  of 

2. _(f)|HQ  y  (  |  Hq,  )  to  zero  and  simplifying  leaves  the  implicit 

equation  for  y^: 


*  & 

£  <  ri9>  -  N  wb  )  -  M(  Ub  +  ab  )  ) 


m  =  1 


t  1  -  -1-j  V  (  r's)  -  N  Jb  )2  ]  -  0  , 

MNa  7  *-> .  m  D 

b  m  =  1 


(  4.34  ) 


where  ic^(a)  =  ^(Q  +  a)  -  ’4*  (a)  which  is  an  increasingly  weak 
A  2 

function  of  a  as  CNR  increases.  Again,  if  was  not  a  function  of  Pb 

through  CNR^,  the  last  term  would  not  appear  and  the  equation  would  reduce 
to  y^  =  the  frame  sample  mean,  as  our  intuition  would  lead  us  to  guess. 
Even  with  the  y^-  relationship,  the  last  term  is  small  for  large  M. 

So,  we  can  take  pb  =  as  an  even  better  approximation  for  the 
background  than  for  a  target  subframe. 


With  y b  and  y  set  to  their  sample  means  for  a  frame  and  subframe, 
respectively,  we  can  complete  deriving  the  processor.  The  conditional 
densities  are  Gaussian  under  Hq  and  ,  so  the  log  likelihood  ratio  is: 


log  A 


:  2 

1  i  b 

2  log  ~*~T 

0  tm 


<  4S>  -  >2 


(  4.35  ) 


We  will  drop  the  term  involving  the  logarithm  of  the  ratio  of  the 
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v 

v 


variances  since  this  term  is  proportional  to  the  logarithm  of  the  contrast 
ratio  and  the  other  term  is  proportional  to  the  contrast  ratio  directly. 
Then  an  equivalent  processor  for  a  Neyman-Pearson  test  is: 

108  '  I  riS>  "  N  I  •  (4’36> 

This  processor  is  very  intuitively  appealing.  We  choose  as  the  target 
subframe  that  subframe  which  exhibits  the  largest  contrast  (positive  or 
negative)  relative  to  the  computed  background  mean.  The  test  is  simply  to 
compare  this  contrast  to  a  threshold. 

Notice  that  this  processor  is  not  quite  in  the  form 

log  =  max{  log  1  (  4.37  ) 

8  m  8 

/*s 

as  we  would  like  since  the  quantity  in  the  braces  contains  U ^  which  is  a 

(g) 

function  of  for  all  m.  However,  since  MN  is  generally  large,  we 

know  the  variance  of  y^  must  be  fairly  small  since  u ^  is  approximately  the 
sample  mean  and  the  variance  of  the  sample  mean  falls  as  1  /  MN.  Because 
of  this,  we  will  take  y  ^  as  exactly  u^,  that  is,  treat  CNR^  as  a  known 
quantity.  Then  the  log  likelihood  ratio  is  in  the  desired  form.  Numerical 

A 

calculations  withy^  =y^  +  3a^  give  indistinguishable  results  with  MN  = 
128  x  128  pixels. 

Although  we  have  justified  assuming  conditional  Gaussian  statistics  for 

C  s ) 

R^  via  a  Central  Limit  Theorem  argument,  we  know  for  CNR  ®  the 

(s) 

random  variable  R  /  N  y  .  has  conditional  chi-squared  statistics  under 
mi 
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rL  with  mean  1  and  2N  degrees  of  freedom.  Let  s  compare  this  chi- 

2  2 

squared  density  to  a  Gaussian  density  with  mean  1  and  variance  a  /Ku  . 

We  will  do  this  comparison  by  calculating  Pp  versus  X  '  for  both  models 

with  X  '=  X  /  N  u  for  Pp  down  to  about  10  .  This  test  allows  us  to 

compare  the  tails  of  the  densities,  which  are  of  most  interest  to  us, 

where  the  Central  Limit  Theorem  arguments  don't  yield  very  accurate 

results.  Figure  4.2  shows  this  comparison  for  M  =  1000.  Although  the  two 

models  converge  more  as  N  gets  larger,  even  at  N  =  100,  the  two  are  still 

-4 

very  far  apart  for  even  fairly  high  false  alarm  rates  like  10  .  Of 

course,  the  problem  is  worse  at  more  reasonable  Pp  values  like  10  ^  which 
are  farther  into  the  tails  of  the  density. 

Although  we  obviously  cannot  use  the  Gaussian  statistics  to  analyze  the 
processor's  performance,  we  will  continue  to  use  the  processor  derived 
from  the  Gaussian  statistics  assumption  because  it  is  so  intuitively 
appealing.  For  analysis  purposes,  we  will  use  an  approximation  derived 
from  the  Chernoff  bounds  [4],  For  CNR  ->  »  ,  numerical  calculations  show 
this  approximation  differs  from  the  exact  results  (the  chi-squared 
statistics)  by  only  a  percent  or  two  [18], 

H.  Range  Only  Processor 

In  this  section  we  will  examine  a  processor  which  uses  all  the  range 
information  available  and  discards  all  the  intensity  information.  For 
each  pixel  we  measure  the  peak  over  all  Q  bins  and  record  the  bin  number 
where  the  peak  occurred,  Q^.  We  ignore  the  intensity  of  this  peak. 

We  derive  the  range  only  processor  using  the  same  techniques  we  derived  in 
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section  C  above  using  the  marginal  coarse  range  probability  mass  functions 
we  derived  in  the  previous  chapter.  As  usual,  we  take  each  pixel  and 
subframe  to  be  independent  of  the  others.  The  resulting  generalized 
likelihood  ratio  test  is: 


A<f> 


max  < 
m  I 
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where  the  denominator  is  evaluated  with  Pr(  )  which  we  in  turn 

derived  from: 

M  N  (  4.39  ) 

p  ?aX  *  El  FI  Pr(  in  ■  "mn  I  H0>  Pr<  :,0b  >  >  >' 

Pr(  “Ob5  m  .  1  n  -  1 


Taking  the  logarithm,  substituting  the  probabilty  mass  function  in, 
differentiating  with  respect  to  Pr(  Xq^  ),  and  setting  the  derivative  to 
zero,  we  solve  for: 


Pr( 


‘q  Q,. 

mn  bmn 


(  4.40  ) 


The  estimate  of  Pr(  ti  ^  )  is  just  the  relative  frequency  of  occurence  of 
anomaly  events  for  the  entire  frame.  If  there  is  no  target  in  the  frame, 
the  estimate  corresponds  exactly  to  the  intuitive  notion  of  the 
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probability  of  no  anomaly  and  the  Law  of  Large  Numbers  [16].  If  there  is 
a  target  in  the  frame,  the  assumption  we  made  of  M  »  1  makes  the  target 
contributions  to  the  estimate  of  the  background  anomaly  probability 
dilute,  and  hence  we  may  use  the  same  intuitive  interpretation  of  the 
equation.  Since  MN  is  the  number  of  pixels  in  a  frame,  it  is  usually 
quite  large,  so  the  estimate  can  be  very  accurate,  particularly  at  small 
background  CNRs.  At  higher  CNRs,  the  estimate  will  be  less  accurate 
(because  the  summation  will  be  smaller  so  statistical  variations  introduce 
larger  relative  errors).  In  the  ensuing  discussion,  we  will  take  the 
estimate  to  be  exact  and  so  treat  the  background  CNR  as  known,  as  we  did 
in  the  intensity  only  receiver. 

Performing  an  identical  maximization  in  the  numerator  gives  the  estimate 

N 

Pr<  a0t)  -T  •  “here  4  5  .  1  E  Sq  Q  '  •  (  4’41  ) 

l<Qt<QT  „  .  1  c 

/\ 

Again,  Pr(  atQt  )  corresponds  to  our  intuitive  notion  of  the  probability 
of  no  anomaly.  Here  is  the  number  of  hits  in  the  range  bin,  from  among 
the  Qrp  potential  target  range  bins,  with  the  most  hits  for  the  m-th 
subframe.  However,  since  N  is  generally  not  very  large  (  <100  ),  the 
estimate  may  not  be  very  accurate.  Hence,  we  cannot  take  this  estimate  to 
be  exact.  Physically,  j  is  the  same  as  it  was  for  the  joint  range  and 
intensity  processor  we  already  examined  and  pictured  in  figure  4.1. 

Without  the  intensity  information  however,  we  choose  Q  in  a  different 
manner  in  equation  4.41.  Now  we  simply  choose  as  Q  the  bin  with  the  most 
hits  from  among  the  potential  target  range  bins  in  the  subframe  picture 
(the  last  line)  in  figure  4.1.  We  can  see  intuitively  now  why  the  range 
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only  processor  must  under-perform  the  joint  range  and  intensity 
processor.  The  joint  processor  takes  into  account  not  only  the  number  of 
hits  per  bin,  but  also  the  average  intensity  for  those  hits.  Since 
anomalies  occur  in  noise  bins,  they  have  intensities  with  mean  1,  which  is 
always  less  than  the  target  bin  with  mean  CNRt  +  1.  Thus,  the  joint 
processor  can  weight  the  range  bin  measurements  to  make  implicit  guesses 
at  whether  or  not  some  hits  correspond  to  anomalies.  Without  the 
intensity  information  to  help  weight  the  decisions,  the  range  only 
processor  will  make  poorer  decisions,  particularly  at  low  CNRs  where 
anomalies  occur  more  frequently. 


If  we  substitute  these  quantities  into  equation  4.38,  take  the  logarithm, 
and  define  the  number  of  hits  at  the  known  background  ranges  as: 


k 

m 


(  4.42  ) 
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(  Q  -  1  )  (  1  -  Pr(  ct0t  )  ) 
(  QT-  1  )  (  1  -  Pr(  a0b  )  ) 
(  Q  -  1  )  Pr( 

1  -  Pr(  aQb  ) 


Notice  that  P^  and  P2  are  functions  of  Pr(  ctQt  )  and  so  are  functions  of 
Because  of  this,  we  can  plot  P^  and  P2  as  functions  of  for  various 
values  of  CNRb  and  get  a  feeling  for  how  sensitive  the  processor  is  to  the 
values  of  P^  and  Figures  4.3  and  4.4  show  such  plots  for  Q  =  100  and 

Qg  <<  Q  ^  Q,p,  a  realistic  condition  for  either  a  large  range  uncertainty, 
L^,  or  a  radar  with  good  range  resolution.  Although  both  parameters  can 
vary  over  a  wide  range,  the  probable  values  for  j  for  reasonable  target 
CNR  values  (>  10  dB)  will  place  the  parameters  in  a  relatively  narrow 
range.  This  being  the  case,  an  approximately  equivalent  processor  for  a 
Neyman-Pearson  test  is  the  much  simpler  formula: 


log  =  j  -  k 

0  u  gm  Jm  m 


(  4.44  ) 


In  fact,  this  processor  is  exactly  optimal  for  known  target  CNR  and  zero 
contrast.  We  will  use  this  processor  for  all  CNRs  since  it  is  a  fair 
approximation  to  the  optimal  processor  and  it  is  possible  to  analyze  its 
performance  without  too  much  difficulty. 


*>1 


Physically,  this  processor  determines  whether  the  presumed  target  range 
bin  or  the  background  range  bins  had  more  hits  and  how  many  more.  As  CNRt 


increases  in  the  subframe  containing  the  target,  log  a 


grows  more 


positive.  As  CNRb  increases  in  a  background  subframe,  the  statistic  grows 
more  negative.  The  processor  looks  for  range  measurements  to  clump  or 
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Figure  4.3 

Parameter  for  Range  Only  Processor 
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Figure  4.4 

Parameter  P-,  for  Range  Only  Processor 
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aggregate  at  the  background  ranges  (the  km  term)  and  somewhere  in  the  Q^, 
potential  target  range  bins  (the  j  term).  It  then  makes  a  decision  based 
on  which  clump  is  bigger  and  by  how  much.  Colloquially,  we  can  say  that 
whereas  the  intensity  only  processor  looks  for  the  maximum  intensity 
contrast,  the  range  only  processor  looks  for  the  maximum  range  contrast. 

I.  Other  Processors 

We  have  derived  processors  for  intensity,  range,  and  joint  intensity  and 
range  measurements.  We  could  also  imagine  incorporating  Doppler 
measurements  into  these  processors.  However,  our  analysis  in  chapter  3 
indicates  we  need  do  no  further  work  on  the  processors  because  of  the 
range-Doppler  duality  we  encountered.  Anywhere  we  treated  range 
measurements,  we  could  substitute  Doppler  measurements,  or  include  joint 
range  and  Doppler  measurements.  The  processors  don’t  change  form, 
although  some  parameters  may  change  values.  Because  of  this  duality,  we 
won't  treat  Doppler  processors  explicitly  until  chapter  6  where  we  will 
consider  how  the  performance  scales  with  changes  in  the  system  parameters 
that  differ  between  range  and  Doppler  radars. 

Now  that  we  have  derived  the  processors  we  will  use,  we  will  derive 
performance  measures  for  them  in  the  next  chapter. 
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V.  Receiver  Performance  Derivations 

Recall  that  we  wish  to  find  Pp  and  P^  for  intensity  only  and  range  only 

receivers  when  they  evaluate  an  entire  frame  of  data.  However,  we  already 

know  we  can  easily  get  the  frame  level  performance  directly  from  the 

subframe  level  performance.  Consequently,  this  chapter  derives  equations 
(s)  (s) 

for  Pp  and  P^  for  a  receiver  measuring  data  only  over  a  single 
subframe.  The  translation  to  frame  level  performance  is  direct.  (For 
notational  simplicity  we  don’t  use  the  (s)  superscripts  in  this  chapter 
nor  the  then  superfluous  m  subscript.  If  there  is  cause  for  confusion,  we 
will  reintroduce  the  subscripts  and  superscripts.) 


A.  Intensity  Only 

Recall  from  the  last  chapter  the  quasi-optimal  intensity  only  processor 
tests  for  subframe  contrast: 


(  5.1  ) 


n  =  1 


If  N  =  1,  we  can  analyze  the  processor  exactly.  The  analysis  is  fairly 
simple  and,  using  the  exact  statistics  derived  in  chapter  3,  the  Pp  and 
P^,  in  terms  of  the  threshold,  are: 
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K  r\  1  ”“34_C  1  X  )  [i  u 

b  )Q  1  (  1  -  e  b  )  u(  1  -  X  1  )  , 

(  5.2  ) 


4  * 


Reasonable  Pp  performance  (  <10  )  requires  \  '  >  1  and  so 

Pp  -  exp[  -  M  b  (  1  +  X’)],  which  gives 

X  ’  ~  ~  (  1  +  log(  Pp)/(ab  Mb)),  and  (  5.3  ) 

PM  a  1  -  exPt  ~  log(  PF  )  ]  =  1  -  (  PF  b 

b 


These  equations  are  very  similar  to  the  equivalent  equations  for  a 
processor  which  measures  intensity  at  the  known  target  range.  The 
equation  for  the  threshold,  X  * ,  in  terms  of  Pp  in  5.3  is  identical  to 
that  for  the  known  range  processor  with  clutter  [14],  The  equation  for  P^ 
in  terms  of  Pp  in  5.3  is  nearly  the  same  as  its  (clutter  free)  analog  for 
the  known  range  processor.  The  only  difference  is  the  exponent  to  Pp 
which,  for  the  known  range  processor,  always  exhibits  an  effective  CNRb  = 


Notice  that  with  zero  contrast,  the  performance  reverts  to  the  guessing 
curve  (Pp  =  Pp)»  as  we  would  expect  for  a  processor  which  searches  for 
contrast.  Similarly,  there  is  an  obvious  asymmetry  between  performance 
for  positive  and  negative  contrasts  of  the  same  magnitude.  This  occurs 
because  the  density  function  is  so  asymmetric  with  a  mode  close  to  the 
origin  and  long  tails  out  to  positive  infinity.  Finally,  notice  that  the 


processor  performance  is  independent  of  Q.  This  occurs  because  X '  is 
large  enough  to  make  the  Q-dependent  terms  1. 


Although  the  single  pixel  processor  is  not  of  much  practical  importance, 
the  ability  to  make  an  exact  analysis  gives  us  information  which  may  carry 
over  to  multipixel  processors  as  well.  We  should  not  be  surprised  if 
multipixel  processors  don’t  perform  well  at  negative  contrasts  and  the 
performance  is  relatively  insensitive  to  the  value  of  Q. 

When  there  are  multiple  pixels  in  a  subframe,  we  have  already  noted  the 
statistics  are  too  complicated  to  analyze  exactly  and  neither  the  Gaussian 
nor  the  chi-squared  limits  are  accurate  enough  for  general  N  and  CNR 
values.  Instead  we  use  an  approximate  performance  measure  derived  from 
the  Chernoff  bound  [4]: 

N 

Pr {  ^  In  2  I  ^  "  exP[  ~  N(  s  Y t(  s  )  -  Yi(  s  )  ] 

n  =  1  _ (  5.4  ) 

•  ^  exp[  i  U  s2  Yi(  s  )  ]  erfc(  VjNsyTcs)  ), 

A 

Yt(  S  )  =  log (  E[  esI  I  H±  ]  ) ,  i  =  t  or  b,  s  >_  0  . 

The  parameter  s  is  free,  but  the  optimal  choice  is  when  y  (  s  )  =  X  . 

The  first  exponential  is  just  the  Chernoff  bound  and  the  other  two  terms 

are  the  correction  term.  For  Ns  y(  s  )  >  9,  we  can  combine  terms  and 


approximate  the  correction  term  by: 


it'*', 


If  the  test  in  equation  5.4  is  <  A  ,  the  only  change  is  in  the  sign 


of  s  in  the  equations,  so  s  <  0. 


For  the  multipixel  intensity  only  processor  we  find  that: 


11 

*  PrM  £  i„  -  N  I  >  »  I  «o  ) 


n  =  1 


-  Pr{  £  In  >  N  yb  (  1  +  X’  )  |  Hq  }  + 


n  =  1 


Pr  {  Y  in  <  N  yb  (  1  -  X’  )  I  H0  }  , 


n  =  1 


JM-Prl  I  Z  P„-N  “b  I  <  *  I  H1  ) 


n  =  1 


-  Pr{  Y,  ^n  >  N  yb  (  1  "  X’  }  i  H1  1  + 


n  =  1 


Y  K  <  n  n ( 1  +  x’ >  i  hi  1  • 


n  =  1 


(  5.6  ) 


The  conditional  semi-invariant  moment  generating  function,  y i(  s  ),  is: 


Yi(  s  )  =  log (  s  [  B(  Q,  ai  -  s  )  -  B(  Q,  -  s  )  ]  )  ,  (  5.7  ) 


i  =  b  or  t, 
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where  B(  *  )  is  the  beta  function: 

r(  x  )  r(  y  ) 

B(  x,  y  )  a  -  .  (  5.8  ) 

r(  x  +  y  ) 


The  s  parameters  in  terms  of  the  thresholds  are: 

y.(  s  )  =  ub(  1  +  X’  )  ,  (  5.9  ) 

0  <  s  ^  for  the  >  X  tests,  s  0  for  the  <  X  tests. 


Going  from  this  point  to  explicit  numerical  results  requires  much  calculus 
and  algebra,  but  there  are  no  new  theoretical  developments,  so  we  will 
skip  these  steps.  In  the  next  chapter  we  will  present  and  discuss  results 
of  actual  computations. 

B.  Range  Only 

We  will  handle  the  range  only  processor  in  a  completely  different  fashion 


from  the  intensity  only  processor.  From  the  last  chapter,  the  processor 
is,  again  for  a  single  subframe,  one  which  searches  for  clumping  of  the 
range  data  in  the  background  and  target  bins: 


H 


1 


i  ~  k 


< 

H, 


(  5.10  ) 


The  Pp  and  values  are: 


Pr{  j  -  k  >_  X  |  Hq  } 
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il  U1C1A. 

Y  Y  Pr(  J  =  j,  K  =  k  I  H0  )  , 


j  =  j  .  k  =  0 
J  Jmin 


j  .  =  max{  0,  X  }  ,  k  =  max{  0,  min{  N  -  j,  X 

Jmin  1  *  ‘  '  max  J 


-  j  )  )  ; 


PM  =  Pr{  j  -  k  <  X  |  Hj  } 


(  5.11  ) 


II  ~  Ill  cl  A 

=  Y  Z  Pr(  J  =  j'  K  =  k  I  Hj  )  , 

k  =  k  .  j  =  0 
min  J 

k  .  =  max(  0,  -  X  }  ,  j  =  max{  0,  mint  N  -  k,  k  +  X  }  }  . 

mm  *  Jmax 


Under  both  hypotheses  we  assume  in  case  of  a  tie,  j  -  k  =  X  ,  we  choose 
hypothesis  H^.  These  are  simple  sums,  but,  of  course,  we  need  to  find  the 
joint  probability  mass  functions  first. 


To  find  the  joint  probability  mass  function,  start  with  a  different  mass 
function  we  can  find  relatively  easily.  We  will  divide  the  Q  range  bins 
into  Qj  potential  target  bins  and  Qg  potential  background  bins  as  before. 

From  the  geometry  discussed  in  chapter  3,  we  know  these  latter  Qg  bins  are 

always  the  most  distant  range  bins  for  a  particular  subframe  and,  for  a 
particular  pixel,  one  of  them  is  the  known  background  range  bin.  Call  J^, 

1  <_  i  <_  Q.J,  the  number  of  hits  in  the  i-th  potential  target  bin,  K  the 

number  of  hits  in  the  correct  background  range  bins  (precisely  which  bin 
is  "correct"  depends  on  the  pixel  number  within  the  frame),  and  L  the 
number  of  hits  in  the  Qg  -  1  background  bins  not  at  the  correct  background 
range.  Under  hypothesis  Hq  we  can  write  the  joint  mass  function 
Pr(  Jj=jp...,Jg  =jg^,K=k,L=l  j  Hq  )  fairly  easily.  It  is  simply  a 
multinomial: 
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Pr(  -  J'i*  •••»  jqt  ~  Jqt*  K  -  k,  L  -  1  I  Hq  )  - 


N! 


h!  ***  v  k!  1! 


Pr(  «0b  )k  (  QB  -  1  )l 


1  "  Pr^  ’Ob  ^ 


N  -  k 


Q  -  1 


-  (1)  *<  «Ob)k  (V) 

(  N  -  k  -  1  )! 


1  ~  Pr<  aOb  } 


,1 N  -  k 


Q  -  1 


(  5.12  ) 


i  !  *  •  •  i  1 

Jl*  JQf 


for  (  jj.  . Jq^  ,  k,  1  }  non-negative  integers  obeying 

+...+  Jq  =  N  -  k  -  1.  Now  define  J  =  max{  |  i  *  1 ,  . . . ,  Q^,  } .  To 

T  1 


convert  the  mass  function  for  the  J.'s  to  a  mass  function  for  J  we  must 

l 


sum  over  all  combinations  of  the  J/s  subject  to  the  stated  constraints 
and  the  additional  constraint  all  ji  <  j: 


Pr(  J  =  j,  K  =  k,  L  =  1  |  Hn  ) 


(  5.13  ) 


Z  Z  Pr(  Ji  =  J'r  •••»  jqt  =  j’qt’  k  -  k,  l  -  i  |  h0  ) 


Ji  =  0"*  v° 


kj  Pr(  a0b 


)k  (N  !  k)  (  Qb  -  1  )X 


1  Pr(  a  Ob 


-iN-k 


Q  -  1 


(Qt) 

C j  *N-k-l 


(Qt)  y  y 

j  ,N-k-l  5  L>  L> 


(  N  -  k  -  1  )! 


i  !  *  ‘  *  i  I 

J.=°  v°  v  V 


j,  +  . . .  +  jp.  -  N  k  -  1 


(Qt) 


V  XT/ 

In  words,  C,  ,  .  is  the  number  of  ways  to  order  the  elements  of  all  the 
J » 


Q,^,  }  such  that  we  meet  the  constraints  stated 


sets  {  |  i  =  1, . . . , 

above.  Now  sum  over  1  to  get  the  mass  function  we  want: 

Pr(  J  -  j,  l  -  k  |  H0  )  -  Q  Pr(  0„  )k 

T  (n  i k)  < «» - 1  >x 

1  =  0 

(Qt) 

Evaluating  ^  is  computationally  cumbersome,  but  once  we  have  the  C- 
values  the  remainder  of  the  calculation  is  straightforward.  Under 
hypothesis  the  procedure  is  similar  though  not  the  exact  dual.  To 
start,  we  assume  the  maximum,  J,  occurs  in  the  correct  (true  target)  range 
bin.  This  assumption  simplifies  the  derivation  and  we  will  shortly  see 
the  error  it  introduces,  for  N  and  CNRt  values  of  interest,  is  very 
small.  Call  the  event  in  which  the  maximum  J.  occurs  in  the  correct  bin 

l 

3  q  and  the  event  in  which  the  maximum  vL  occurs  in  the  wrong  bin  3  ^ • 

Then  the  mass  function  we  start  from  is: 


1  -  Ob  > 

Q  -  1 


N-j-k 

Pr(  J  =  j,  K  =  k,  30  |  Hj  )  =  Y,  Pr(  J  =  j,  K  =  k,  L  =  1,  SQ  |  H]>  ) 

1  =  0 


(  5.15  ) 


The  difference  between  C’  and  C  is,  since  we  have  8q,  there  is  one  less 
sura  in  the  sequence  of  summations  and  there  are  only  Q^,  -  1  bins  to 
consider.  For  N  =  10  and  CNRt  =  10  -0,  both  fairly  small  values, 


Pr(  Bq  |  )  =  0.99966.  So,  8g  is  essentially  always  true  for  cases 

of  interest  to  us. 

With  these  probability  mass  functions  we  can  compute  Pp  and  P^  directly 
using  the  equations  in  5.11.  However,  Pp  is  not  a  continuous  function  of 

(s) 

the  threshold  X  since  the  test  statistic  log  A  is  discrete.  We 

gm 

touched  on  this  problem  at  the  beginning  of  chapter  4  and  noted  we  could 
make  Pp  a  continuous  function  of  X  by  randomizing  the  decision  rule  with 
a  well  known  procedure.  We  will  use  this  practice  when  computing  Pp  and 
P^  for  the  results  in  the  next  chapter  so  we  can  use  Pp  and  P^  values  in 
the  range  processor  identical  to  those  we  used  in  the  intensity 
processor.  This  will  make  it  much  easier  to  do  comparisons  between  the 
two  processors  without  introducing  corrections  to  make  the  comparisons 
seem  fair. 

Now  we  have  processors  for  range  only  and  intensity  only  receivers  and  we 
have  methods  to  compute  their  performance  statistics.  Actually  performing 
the  computations  can  be  a  complicated  numerical  procedure,  but  detailing 
the  procedure  would  not  contribute  to  the  main  point  of  this  thesis  and 
there  are  no  novel  numerical  techniques  necessary.  Consequently,  we  will 
turn  now  to  simply  examining  the  results  of  these  computations. 


VI.  RECEIVER  PERFORMANCE  RESULTS 


We  have  taken  our  theoretical  development  from  the  single  pixel  speckle 
target  IF  waveform  through  a  multipixel,  multidimensional  probability 
density  function  and  derived  quasi-optimal  processors  for  two  important 
receivers,  the  range  only  and  intensity  only  measurement  receivers.  We 
have  shown  how  to  compute  the  P^-Pp  performance  curves  for  these  receivers 
in  the  last  chapter.  These  theoretical  developments  are  all  we  need  to 
meet  the  original  objectives  of  this  thesis.  Now  we  will  present  some 
actual  numbers  to  demonstrate  the  relationships  embodied  in  the  processor 
and  performance  equations  we  derived  in  the  last  two  chapters.  In  this 
chapter  we  will  first  derive  relationships  for  how  the  CNR  (target  or 
background)  scales  with  changes  in  various  system  parameters  and  then 
actually  compute  the  performance  curves  for  some  cases  of  these  processors 
and  compare  the  range  and  intensity  only  receivers.  As  in  previous 
chapters,  we  will  perform  our  calculations  only  for  a  ranging  radar  and 
comment  afterwards  on  how  to  carry  the  results  over  to  Doppler  and  range- 
Doppler  receivers. 

A.  System  Parameter  Interdependence 

Our  equations  for  system  performance  involve  many  parameters:  M,  the 
number  of  target  sized  subframes  in  one  frame  or  field  of  regard;  N,  the 
number  of  pixels  per  subframe  or  on  target;  Q,  the  number  of  range  bins 
per  pixel;  CNRt,  the  target  CNR;  CNR^,  the  background  CNR;  and 
which  we  will  define  as  the  contrast  ratio,  CNR^^/CNR^.  These  parameters 
are  not  all  independent  and  some  depend  on  other  system  parameters  which 
do  not  appear  explicitly  in  the  performance  equations.  However,  the  radar 
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equation,  3.17,  determines  how  these  parameters  interact  to  affect 
performance: 


n  P«j> 

hv  B 
o 


-2a  L. 


(  6.1  ) 


For  simplicity  of  presentation,  we  will  use  a  =  0  and  e0pt  =  1  and  keep 
fixed  n/hv^  and  target  and  background  reflectivity  so  that: 


Now  consider  how  the  choice  of  the  parameters  listed  in  the  first  sentence 
above  affect  the  CNR  through  the  radar  equation,  6.2,  and  vice-versa.  To 
do  this,  we  will  first  review  the  models  for  the  target  and  its 
environment  and  decide  which  elements  of  these  models  we  will  have 
reasonable  control  over  and  which  we  won't. 

We  will  assume  a  fixed  search  volume  and  target  size.  The  fixed  search 
volume  corresponds  to  a  fixed  angular  uncertainty,  n  steradians,  and 
fixed  range  uncertainty,  Ly  meters.  These  assumptions  simply  say  the 
radar  designer  has  no  control  over  the  size  of  the  target  or  where  it 
resides  in  space.  For  any  given  problem,  however,  the  radar  designer 
knows  the  fixed  values  of  ft  ,  Lu,  and  the  target  size,  ft  steradians. 
Notice  we  jumped  from  knowing  the  target  size  to  knowing  the  target 
angular  extent  ft  .  Formally,  this  requires  knowing  the  target  range,  L, 
which  we  don’t  know.  However,  if  we  know  there  is  only  a  small  relative 
difference  between  the  minimum  and  maximum  L  values,  will  be 
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approximately  constant  for  any  L  in  the  range  of  permitted  values.  We 
have  assumed  a  known  n  all  along  by  assuming  a  known  N.  A  system  where 
is  not  known  would  require  treating  N  as  an  unknown  parameter  in  the 
fashion  of  chapter  3. 


For  a  given  Lu,  we  have 


(  6.3  ) 


where  1  is  the  bin  size  or  range  resolution  in  meters  and  tf  is  the  bin 
size  or  radar  pulse  width  in  seconds.  Since  all  the  performance  equations 
depend  only  on  Q,  Ly  does  not  enter  the  equations  explicitly.  For  a  given 
q.,  the  value  of  N  depends  on  the  angular  resolution,  Aft  steradians: 


(  6.4  ) 


however,  M  remains  constant  at 


(  6.5  ) 


Changing  the  radar  parameters  has  no  effect  on  M.  Because  this  is  true, 
and  because  performance  depends  on  M  in  such  a  predictable  fashion,  we 
will  not  concern  ourselves  too  much  with  M.  We  will  fix  M  =  1000,  as  we 
did  in  chapter  5  for  a  single  pixel,  for  the  remainder  of  our  work.  This 
corresponds  to  a  frame  of  about  30  by  30  subframes. 


We  should  mention  there  is  one  complication  with  this  assumption.  Some  of 


■‘.r.r.  r,  •- 


v-  «-.v. wi'.-.v  - J 


our  derivations  assumed  a  large  M  value.  Although  we  didn't  define  what 
we  meant  by  large,  cases  of  M  10  are  probably  too  small  to  allow  the 
approximations  to  remain  valid.  We  won't  deal  with  this  problem,  but  we 
should  at  least  recognize  it  exists.  Notice  that  N  depends  on  Aft  which 
in  turn  depends  on  the  receiver  aperture  Ap,  whereas  Q  depends  on  tf  which 
in  turn  sets  the  matched  filter  bandwidth  B. 

Now  let's  examine  the  B  or  Q  dependence  more  carefully.  We  will  let  B  = 
l/tr  since  we  used  a  matched  IF  filter.  Pp  is  the  transmitter  peak  power 
and  it  is  related  to  the  transmitter  average  power,  P^,  by 

T  PT 

PT  “  PA  -  -  =  PT  fcr  =  PA  T‘  (  6’6  } 

fcr  ® 

where  T  is  the  pixel  dwell  time  or  inter-pulse  time,  and  1/T  is  the  pulse 
repetition  frequency  (prf).  There  is  a  maximum  pulse  repetition  frequency 
in  a  ranging  radar  that  is  set  by  the  desire  for  unambiguous  range 
information  [1],  [7],  If  we  want  unambiguous  range  detection  out  to  range 
L,  we  must  have  T  >  2L/c;  for  L  =  3  km  this  translates  to  T  >  20  usee. 

We  will  assume  T  takes  on  its  minimum  value,  since  larger  values  lead  to 
longer  frame  imaging  times,  MNT. 

If  we  assume  our  laser  is  average  power  limited,  P^  is  constant.  Since  T 
is  also  constant,  Pp  /  B  is  a  constant  by  equation  6.6.  Hence,  changing  Q 
does  not  affect  the  system  CNR. 

If  we  assume  our  laser  is  peak  power  limited,  Pp  is  constant  and 
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(  6.7  ) 


so  the  system  CNR  is  proportional  to  1/Q.  In  this  case  it  is  useful  to 
define  a  quantity  as  the  largest  number  of  range  bins. 


C  max ) 

corresponding  to  maximum  bandwidth,  B  ,  operation  of  the  radar.  In 


terms  of  these  maximum  parameters  we  have: 


B  =  -1 

2L  Q^max) 

u  x 


B 


(max) 


.(max) 


(  6.8  ) 


.(max) 


The  value  Q  exists  because  we  have  assumed  all  targets  and  background 

.(max) 


are  not  range  resolved  within  a  pixel.  If  we  make  tr  too  small  (Qv 

too  large),  we  will  start  to  resolve  the  fine  range  structure  of  the 

.(max) 


target.  There  is  also  an  indirect  N  dependence  in  Qv  .  For  a  given 
target,  if  N  is  small,  each  pixel  will  cover  a  larger  area  of  the  target 
which  leaves  more  opportunities  for  the  target  to  be  range  resolved.  For 


large  N,  we  have  a  larger  Q^max\  We  won't  deal  with  this  weak  N 


dependence  in  this  thesis. 


Turning  to  the  N  parameter,  we  use  Fourier  diffraction  theory  [21]  and 
equation  6.4  to  obtain 


1.2 


1.2 


Afl 


(  6.9  ) 


where  the  factor  1.2  arises  from  the  Airy  disk  radius,  tt  ,  factors  of  2, 


(max) 


etc.  If  we  define  a  maximum  aperture  ApV  '  and  a  corresponding  maximum 
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1.2  X2  N(max)  N(B,a’t) 
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_ _ _  _  ^(max)  _ _ _ 

N(max)  “  AP  M(max) 


Finally,  in  terms  of  the  CNR  at  maximum  aperture  and  maximum  bandwidth, 
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we  have  the  scaling  laws 
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These  relationships  are  important  for  making  reasonable  comparisons 
between  two  processors  with  different  resolution  parameters  imaging  the 
same  target  and  environment.  We  will  return  to  these  equations  when  we 
compare  processors  with  different  resolutions  and  in  the  last  section  of 
this  chapter.  For  individual  processors  with  fixed  N  and  Q  values, 
however,  there  is  no  need  to  bother  with  the  scaling  factors. 


In  all  the  following  calculations,  we  will  use  N^max)  =  40  pixels  and 
q(max)  ^ 


10000  bins. 


B.  Intensity  Processor 

Recall  the  quasi-optimal  intensity  only  processor  searches  the  frame  for 
the  subframe  which  exhibits  the  most  contrast  relative  to  the  computed 
background  intensity.  If  that  contrast  exceeds  a  threshold,  the  processor 
declares  a  target  present  in  the  subframe  with  the  maximum  contrast: 


i  -  n  u.  |  }  ;  a 

mn  b  1  < 


(  6.13  ) 


n  =  1 


Using  the  equations  we  derived  in  chapter  5,  we  can  plot  versus  CNRt 

for  this  processor  at  given  Pp,  N,  M,  Q,  and  values.  Figure  6.1 

shows  the  frame  level  performance  results  for  N  =  10  and  20  pixels,  M  = 

1000,  and  Q  =  10.  For  each  plot  there  are  two  sets  of  two  curves  each. 

The  top  set  is  for  =  +5  dB  and  the  bottom  set  is  for  =  +10 

“6 

dB.  For  each  set,  the  top  curve  (dotted  line)  is  for  Pp  =  10  and  the 

_3 

bottom  curve  is  for  Pp  =  10  .  For  ®  dB,  all  the  curves  are 

indistinguishable  from  P^  =  1.0.  The  value  of  Q  =  10  is  fairly  arbitrary 
because,  as  we  shall  see  shortly,  Q  is  not  very  important  to  the 
performance  of  an  intensity  only  processor.  We  expected  this  Q- 
insensitivity  from  our  single  pixel  work  in  chapter  4. 


The  method  used  to  generate  the  numbers  in  these  curves  has  introduced  two 
quirks  in  these  plots  which  the  reader  will  easily  notice.  In  the  N  =  10 
plot,  figure  6.1  (a),  the  lines  for  Pp  *  10  ^  don't  extend  to  CNRfc  =  30 
dB.  This  is  not  a  major  problem  since  the  curve  has  approximately  reached 
its  asymptotic  value  by  the  time  the  data  stops  anyway.  In  figure  6.1  (b) 
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for  Pp  =  10  and  =  5  dB,  there  is  a  bump  in  the  curve  at  very  low 

CNRt  values.  This  arises  because  the  Chernoff  bounds  don't  work  well  at 
estimating  large  probabilities.  This  is  not  a  real  problem  since  this 
region  corresponds  to  very  poor  performance  anyway. 


The  first  thing  to  notice  on  these  plots  is  the  asymptotic  behavior  for 

large  CNR^.  As  CNRt  — >  00  ,  the  does  not  go  to  zero.  If  we  use  the 

unit  mean  chi-squared  random  variable  with  2N  degrees  of  freedom  to  model 

R^(s)/g[R^(s) |h^]  for  very  large  CNR^  values,  i  =  t,  b,  we  can  more 

easily  understand  this  behavior.  Let's  define  a  random  variable  I  ^ : 

J  mn 


l(°)  _ 


I 

mn 


- - W 

E[  I 


mn 


(  6.14  ) 


where  Imn^  is  a  conditionally  chi-squared  random  variable  with  2  degrees 
of  freedom  and  mean  1  under  Hq,  mean  under  H^.  The  sum  of  N 

conditionally  chi-squared  random  variables  with  2  degrees  of  freedom  is 
another  conditional  chi-squared  random  variable  with  2N  degrees  of  freedom 
[16].  Now  let  us  compute  the  relative  distance  between  the  two 
conditional  means: 


(  6.15  ) 


Notice  for  a  fixed  contrast  ratio,  the  distance  is  fixed.  As  we  increase 


,  the  densities  converge  to  these  chi-squared 


CNR^  with  constant 
densities  and  the  distance  between  the  densities  for  the  intensity 
processor’s  decision  statistic  under  Hq  and  don't  spread  farther 
apart.  Consequently,  the  two  hypotheses  don't  get  any  easier  to 
distinguish  as  CNR  increases  at  constant  contrast,  so,  approaches  a 
nonzero  asymptotic  value. 

Physically,  we  are  at  speckle  limited  performance.  The  CNRs  are  so  large 
there  is  virtually  no  LO  shot  noise  to  contend  with  and  the  speckle 
induced  intensity  fluctuations  alone  make  the  target  and  background 
difficult  to  distinguish. 


Equation  6.15  also  helps  explain  what  happens  to  performance  as  a  function 
of  £cnr.  For  £cnr  =  1,  no  contrast,  the  distance  is  zero.  Since, 
intuitively,  our  processor  searches  for  contrast,  it  is  logical  it  would 
fail  completely  when  there  is  no  contrast.  However,  equation  6.15  also 
helps  explain  the  asymmetry  between  performance  for  positive  and  negative 
contrasts  of  the  same  magnitude.  For  £  >  1,  d  grows  linearly  with 

^  CNR  *  However,  for  £qjR  <  1,  d  ->  1  as  contrast  increases  in 
magnitude.  This  limiting  value,  a  result  of  the  asymmetry  in  the 
densities  around  their  means,  limits  the  performance  for  even  very  large 


negative  contrasts. 


Physically,  for  negative  contrasts  the  mean  target  intensity  is  lower  than 
the  mean  background  intensity.  However,  the  speckle  induced  intensity 
fluctuations  in  the  background  measurements  often  give  intensities  less 
than  the  mean  and  close  to  the  target  intensities.  This  makes  the  two 
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virtually  impossible  to  distinguish. 


Now  let  us  examine  the  effect  of  Q  and  N,  i.e.,  range  and  spatial 
resolution,  on  the  performance.  We  will  plot  versus  Q  or  N  for  fixed 
Pp,  M,  and  for  both  laser  power  models.  Figure  6.2  plots  P^  versus 

Q  for  N  =  10  and  N  =  20  at  a  CNR  of  16  and  20  dB.  Figure  6.2(a)  is 

for  the  constant  average  laser  power  model  and  figure  6.2(b)  is  for  the 
constant  peak  laser  power  model.  Figure  6.3  plots  P^  versus  N  for  Q  =  10 
and  1000  at  CNR  0f  15  an(j  20  dB.  Figure  6.3(a)  is  for  the  constant 

average  laser  power  model  and  figure  6.3(b)  is  for  the  constant  peak  laser 
power  model.  All  plots  hold  constant  M  =  1000,  Pp  =  10  and  Cqjr  *  +10 
dB.  These  figures  incorporate  the  performance  corrections  for  different  N 
and  Q  values  we  derived  in  section  A  of  this  chapter. 


M 


For  the  constant  average  laser  power  model,  the  performance  is  almost 
completely  independent  of  Q  as  we  expected  from  the  single  pixel 
equations.  To  the  extent  there  is  some  dependence,  the  performance  gets 
poorer  as  Q  increases  because  of  the  increased  probability  of  anomaly. 

For  the  constant  peak  laser  power  model,  there  is  a  Q  dependence,  but  it 
is  due  to  the  CNR  dependence  on  Q.  As  we  decrease  Q  from  the  full 
bandwidth  value  Q^max^|  the  CNR  rises  and  performance  improves  until  it 
saturates  at  the  speckle  limited  values.  The  slight  drop  in  performance 
at  very  low  Q  values  for  N  =  10  is  not  real.  It  is  a  result  of  the  lack 
of  data  at  high  CNR  values  mentioned  in  describing  figure  6.1(a).  (When  I 
ran  out  of  data  for  higher  CNRfc  values,  I  assumed  the  final  data  point  had 
reached  its  asymptote.  This  assumption  is  not  quite  valid,  as  is  obvious 
from  figure  6.1(a).) 
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From  the  form  of  the  Chernoff  bound,  we  expect  the  performance  to  get 
exponentially  better  with  increasing  N  and  figure  6.3  bears  this  out  for 
both  models.  Notice  in  the  constant  average  laser  power  model  the  4  dB 
change  in  CNR  makes  a  bigger  performance  difference  than  the  2  order 
of  magnitude  change  in  Q.  For  the  constant  peak  laser  power  model,  the  Q 
=  10  curves  coincide  for  both  CNR  values  because  there  is  a  30  dB 
increase  in  CNR  from  Q^max^  =  10000  to  Q  =  10,  so  the  processors  are  in 
their  speckle  limited  regimes  regardless  of  the  4  dB  difference  in  the 
CNR^O^  values.  For  the  Q  =  1000  case,  there  is  only  a  10  dB  increase  in 
CNR^  over  the  full  bandwidth  CNR,  so  the  4  dB  change  in  CNR  makes  a  small 
difference.  These  curves  are  still  very  close  to  each  other  and  the 
speckle  limited  performance  regime. 

C.  Range  Processor 

Recall  the  quasi-optimal  range  only  processor  determines  the  difference 
between  the  number  of  hits  in  the  presumed  target  range  bin  and  the  number 
of  hits  in  the  known  background  range  bins  for  each  subframe.  If  the 
maximum  of  these  differences  exceeds  a  threshold,  the  processor  declares  a 
target  present  in  the  subframe  where  the  maximum  occurred: 

H1 

max{  1  -  k  }  *  \  .  (  6.16  ) 

m  m  s 
m  „ 

H0 

Using  the  equations  we  derived  in  chapter  5,  we  can  plot  versus  CNRt, 

Q,  and  N  for  the  range  only  processor  just  as  we  did  above  for  the 
intensity  only  processor.  (Since  Qg  enters  the  probability  mass  functions 
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for  j  and  k,  we  would  expect  that  performance  is  also  a  function  of  Qg. 

In  all  ensuing  examples  we  assume  a  square  target  with  each  horizontal 

line  of  pixels  on  the  target  at  a  different,  adjacent  range  bin  number 

from  the  adjoining  line  of  pixels.  Under  this  condition,  N  uniquely 

determines  Qg:  Qg  =  fV»l  .)  Figure  6.4  shows  frame  level  versus  CNRt 

performance  for  N  =  10  and  20  pixels,  M  =  1000,  and  Q  =  10.  Each  plot  is 

—3  —6 

for  Pp  values  of  10  and  10  (dotted  lines).  For  each  Pp  value  there 
are  curves  for  =  +10,  0,  and  -10  dB. 


For  this  processor,  P^  -»  0  as  CNRt  00  for  a  fixed  ^  CNR  *  ^oth  CNRt 
and  CNRg  00  ,  the  probabilities  of  anomaly,  the  probability  of  choosing 
the  wrong  range  bin,  for  target  and  background  0,  so  E[  j  -  k  |  Hq  ]  -> 
-N  and  E[  j  -  k  |  ]  ->  +N.  But  var[  j  -  k  |  H.  ]  ->  0  for  i  =  0,  1, 

so,  if  we  define  a  d  parameter  as  in  equation  (6.14),  we  find  d  -»  ®  . 

The  two  hypotheses  are  easy  to  distinguish  and  the  performance  improves 
with  increasing  CNRt  at  constant  contrast.  Physically,  since  the 
probability  of  anomaly  goes  to  zero,  the  range  estimates  are  always 
perfect  and  so  there  is  infinite  range  contrast. 


We  can  also  explain  the  behavior  as  a  function  of  by  examining  the 

processor.  The  processor  looks  for  range  measurements  from  individual 
pixels  to  concentrate  or  clump  together  in  one  of  two  range  bins,  namely 
the  target  or  the  background  bins.  As  CNR^  increases,  the  probability  of 
anomaly,  Pr(  a  ^  ),  drops  so  the  clumping  is  more  pronounced  in  bin 
for  the  target  or  bin  Qg  for  the  background.  If  we  fix  CNRt  and  vary 
’CNR’  we  change  CNRg.  As  Cqjr  increases  for  a  fixed  CNRg,  CNRg 

drops  and  Pr(  a  jg  )  increases,  so  performance  grows  poorer.  This  is 


precisely  what  figure  6.4  shows:  performance  for  a  fixed  CNR^  varies 
inversely  with 


Now  let  us  examine  the  range  only  processor's  performance  as  a  function  of 
Q  and  N.  We  will  plot  versus  Q  and  N  for  Pj,  =  10  M  =  1000,  and 
^ CNR  =  for  ^oth  laser  power  models.  Figure  6.5  plots  P^  versus  Q 

for  N  =*  10  and  20  pixels  at  CNR^^  =  16  dB  and  20  dB  (dotted  curves). 
Figure  6.5(a)  is  for  the  constant  average  laser  power  model  and  figure 


6.5(b)  is  for  the  constant  peak  laser  power  model.  Figure  6.6  plots  P^ 
versus  N  for  Q  =  10  and  1000  at  CNR^^  =  16  dB  and  20  dB  (dotted  curves) 
under  the  constant  average  laser  power  model. 


Figure  6.5(a)  shows  performance  for  the  constant  average  laser  power  model 
at  first  increases  as  Q  increases  and  then  falls  slightly.  There  are  two 
effects  to  consider  to  explain  this  behavior.  As  Q  increases,  so  does 
Q^,.  Since  anomalies  occur  in  random  bin  numbers,  the  larger  is  Q^,,  the 
less  likely  it  is  that  a  random  clumping  of  anomalies  will  be  large  enough 
to  mistake  for  a  target.  This  causes  performance  to  improve.  On  the 
other  hand,  as  Q^,  increases,  Pr(  a^t  )  increases  approximately 
logarithmically  [4]: 

Pr(  alt  )  “  at  (  l0g(  Q  >  -  TQ  +  °*577  )  •  (  6.17  ) 

which  in  turn  reduces  the  processor's  performance.  These  two  effects 
combine  to  produce  the  behavior  in  figure  6.5(a).  Under  the  constant  peak 
laser  power  model,  the  increase  in  CNRt  as  Q  falls  overwhelms  the  weaker  Q 
dependence  from  the  statistics.  So,  performance  gets  better  very  quickly 
for  decreasing  Q. 


Figure  6.6  shows  that  the  range  only  processor's  performance,  under  the 
constant  average  laser  power  model,  as  a  function  of  N  improves 
approximately  exponentially  in  N.  Since  we  are  dealing  with  sums  of 
random  variables  (we  sum  the  number  of  hits  in  each  bin  over  the 
subframe),  the  Chernoff  bound  applies  to  the  tails  of  the  range  densities 
(although  it  would  be  difficult  to  generate  one  explicitly  for  this 
density).  The  Chernoff  bound  exhibits  this  exponential  dependence  on  N 
for  any  density.  Like  the  intensity  only  processor,  we  see  the  effect  of 
the  change  in  Q  of  two  orders  of  magnitude  is  less  important  than  a  change 
of  4  dB  in  CNRt(0).  There  is  no  curve  for  the  constant  peak  laser  power 
model  because,  for  these  parameter  values,  the  10  dB  increase  in  CNR^  due 
to  reducing  Q  from  10000  to  1000  is  enough  to  drive  all  the  curves  down 
below  10  6. 

D.  Performance  Comparisons 

Now  we  want  to  compare  detection  performance  of  the  range  only  and 
intensity  only  processors.  There  are  many  ways  to  perform  these 
comparisons,  some  probably  more  suitable  than  others  for  any  given 
problem.  I  have  selected  two  methods  of  presentation.  Both  methods  keep 
in  mind  our  goal  of  expressing  system  performance  in  terms  of  system 
parameters  and  resolving  power  in  the  different  measurement  dimensions. 

The  first  method  is  closely  related  to  the  way  we  presented  results  for 
the  processors  individually  in  sections  B  and  C  above.  We  will  plot  the 
quantity  log^  '''nten3:'"t^  /  p^(^an8e)  )  as  a  function  of  Q  and  N, 
where  the  superscripts  indicate  the  processor  from  which  we  derived  the  P^ 
value.  Figure  6.7  shows  such  a  plot  for  the  constant  average  laser  power 


=  16dB  and  18 


model  with  M  =  1000,  Pp  *  10~6,  ^CNR  =  +10  dB,  and  CNRt(0) 
dB.  A  positive  value  indicates  the  range  processor  outperforms  the 
intensity  processor.  It  may  be  difficult  to  see  on  these  plots,  but  the 
intensity  processor  does  slightly  out-perform  the  range  processor  for  low 
N  and  Q  values.  We  saw  both  processors  improve  exponentially  with 
increasing  N,  but  these  plots  indicate  the  range  only  processor  improves 
faster  than  the  intensity  only  processor,  hence  the  surface  always  slopes 
up  for  increasing  N.  Physically,  this  corresponds  to  performance 
increasing  more  quickly  as  a  result  of  range  clumping  in  the  range  only 
processor  than  as  a  result  of  averaging  out  the  speckle  intensity 
variations  in  the  intensity  only  processor.  The  concavity  of  a  plot  for  a 
cut  along  a  constant  N  value  reflects  the  behavior  in  figures  6.5(a)  and 
6.2(a).  Since  the  intensity  only  processor  is  only  very  weakly  dependent 
on  Q,  the  result  is  mostly  due  to  the  change  in  performance  for  the  range 
only  processor  as  a  function  of  Q.  We  have  already  discussed  this 
behavior  in  section  C. 

This  method  shows  quite  graphically  how  the  two  processors  perform 
relative  to  each  other,  however,  we  need  a  different  plot  for  each  CNR 
and  we  get  no  indication  of  the  absolute  P^  for  any  point.  Also,  the 
strong  dependence  of  P^  on  Q  under  the  constant  peak  laser  power  model 
makes  it  difficult  to  construct  similar  plots  for  that  model. 

Consider  the  other  display  method.  Suppose  we  set  certain  performance 
requirements:  require  meeting  a  Pp  and  P^  criteria  at  a  given  CNRt  and 
’CNR*  Then  we  can  simply  check  whether  or  not  a  processor  meets  these 
requirements  as  a  function  of  Q  and  N.  Figures  6.8  and  6.9  show  several 
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Figure  6.S 

Range  Only- Intensity  Only  Processor  Performance  Trades 
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Figure  6 .9 

Range  Only-Intensity  Only  Processor  Performance  Trades: 

Constant  P^  Model 

(a)  High  Performance,  High  Contrast 

(b)  Low  Performance,  Low  Contrast 


displays  for  M  =*  1000  and  various  performance  requirements  under  the 
constant  average  and  peak  laser  power  model,  respectively.  It  is  quite 
easy  to  generate  such  displays  for  a  given  set  of  performance 
requirements.  An  "X”  indicates  the  range  processor  meets  the  performance 
criteria  for  that  Q  and  N  while  an  "0"  indicates  the  intensity  processor 
meets  the  criteria.  From  the  display  a  system  engineer  can,  for  instance, 
determine  what  combinations  of  N  and  Q  will  give  a  range  or  intensity 
processor  which  meets  the  necessary  performance  requirements. 

Figure  6.8,  for  the  constant  average  laser  power  model,  shows  the  more 
significant  impact  of  laser  aperture  size,  Ap,  (reflected  in  the  N  value) 
on  performance  relative  to  the  resolution  parameter  Q.  Generally  the 
processor  either  meets  requirements  for  a  given  N  or  not,  regardless  of 
Q.  Figure  6.8(a)  is  for  relatively  high  performance  requirements:  Pp  = 
99.9%  at  Pp  =  10  ^  and  fairly  low  CNRt^  values.  At  16dB,  we  find  the 
point  N  =  20  pixels,  Q  =  10  bins  where  the  intensity  only  processor  out¬ 
performs  the  range  only  processor.  At  higher  Q  values,  300,  the 
intensity  only  processor  fails,  but  the  range  only  processor  satisfies 
the  requirements  for  Q  _>  30.  If  we  increase  the  CNR^^  by  only  4  dB,  we 
can  meet  requirements  with  either  processor  as  long  as  N  ^  20. 

In  figure  6.8(b)  we  relax  the  performance  requirements  somewhat  to  Pp  = 

95%  at  Pp  =*  10  ,  but  we  reduce  the  contrast  to  only  5  dB.  Here  the 

intensity  only  processors  need  at  least  N  =  30  pixels  to  meet  the  relaxed 
requirements  because  they  perform  so  poorly  at  low  contrasts.  The  range 
only  processors  perform  well  enough,  however,  that  for  a  4  dB  increase  in 
CNRt^\  we  can  cut  angular  resolution  in  half,  from  N  =  20  pixels  to  N  = 


10  pixels  and  still  meet  requirements. 

In  figure  6.9(b),  for  the  constant  peak  laser  power  model,  the  increase  in 
CNR  for  reduced  Q  allows  the  processors  to  meet  performance  at  lower  N 
values.  This  effect  is  strong  enough  to  change  the  Q  dependence  for  the 
range  only  processor  and  make  it  perform  better  for  small  Q  values.  We 
saw  this  phenomenon  in  figure  6.5(b)  as  well. 

E.  Other  Processors 

Now  it  is  worthwhile  to  say  a  few  words  about  other  types  of  processors 
and  how  to  apply  these  results  to  them.  We  derived  dualities  for  Doppler 
and  range-Doppler  processors  in  the  last  two  sections  of  chapter  3.  We 
can  apply  the  results  from  this  chapter  to  Doppler  receivers  as  a  result 
of  those  dualities,  but  we  must  take  care  in  how  we  do  the  comparisons 
because  of  the  system  parameter  interdependencies  we  worked  with  at  the 
beginning  of  this  chapter.  In  general,  we  have  to  consider  each 
transmitter  waveform  separately  and  determine  the  effect  on  scaling 
parameters,  particularly  for  range-Doppler  systems  where  complex  waveforms 
are  often  the  norm. 

For  a  Doppler  only  processor,  we  will  use  the  long  Gaussian  pulse  waveform 
we  used  in  deriving  the  duality  of  chapter  3.  We  must  determine  how  CNR 
scales  with  changing  fr>  frequency  resolution,  and  hence  Q,  for  this 
waveform.  We  still  have  a  pulsed  Gaussian  transmitter  waveform,  so  the 
equations  in  6.6  above  still  hold.  Now,  however,  the  frequency  resolution 
is  approximately  f r  *  1 / t r  and  the  bandwidth  is  approximately  B  =  1/t  = 
f  .  If  F  is  the  frequency  uncertainty  interval  width,  Q  =  F  /f  =  F  /B. 


In  an  exact  analogy  with  the  ranging  radar ,  we  will  assume  the  problem 
statement  fixes  F^.  In  the  ranging  radar,  we  could  make  tr  arbitrarily 
small  without  worrying  about  L^.  In  the  Doppler  radar,  we  must  be  a 
little  concerned  about  making  f  too  small.  Although  the  design  of 
Fourier  transform  devices  like  surface  acoustic  wave  devices  (SAWs)  is 
somewhat  afield  of  this  thesis  topic  [22],  there  is  a  time-bandwidth  (TW) 
product  constraint  on  F^t^..  However,  high  time-bandwidth  product  devices 
are  available  [22]  so  we  will  assume  F^  is  small  enough  that  we  can 
increase  the  frequency  resolution  without  affecting  F  .  Under  this 
condition. 


PT  PT  Q 

— -  =  P.  T  = 

B  Fu 
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Now  we  can  write  how  CNR  depends  on  Q  (the  N  dependence  has  not  changed 
from  the  range  radar  analysis): 


P^,  limited 

P.  limited 
A 


(  6.19  ) 


where  i  *  t  or  b. 


Here  CNR^^  still  corresponds  to  CNR^  at  N  =  N^max^  and  Q  =  Q^max\  Now, 
however,  by  equation  6.18,  Q^max^  corresponds  to  minimum  IF  bandwidth 
instead  of  maximum  IF  bandwidth.  If  we  defined  a  minimum  Q  value,  Q  '  , 

and  defined  CNRi^0)  at  N  -  N^max)  and  Q  =  Q(mln\  then  the  scaling  law  for 
CNRi  would  be  the  same  as  equation  6.12  for  the  ranging  radar  with  Q^max^ 
replaced  by  Q^min\  Then  the  plots  of  figure  6.9  would  still  apply  to  the 


Doppler  only,  constant  model  analysis  except  the  Q  labels  would 
reverse,  running  left  to  right  from  Q  =  10000  to  Q  =  10  (i.e.,  take  a 
mirror  image  of  each  plot  through  an  axis  at  log^CQ)  =  2.5).  However, 
q(max)  gtj^  corresponds  to  maximum  Doppler  resolution  which  is  a  more 
physically  meaningful  quantity  than  Q*'n”‘n\  Analogous  to  the  range  radar 
case,  there  is  a  limit  to  how  fine  we  can  make  the  Doppler  resolution.  If 
we  make  the  resolution  too  fine,  we  will  start  to  resolve  the  fine  Doppler 
structures  of  the  target  like  vibrations  and  rotations  of  various  parts  of 
the  target  within  a  pixel.  This  violates  the  unresolved  Doppler 
assumption  we  used  in  the  derivations.  Consequently,  we  will  use  the 
definition  of  CNR^^  in  terms  of  full  aperture  and  minimum  bandwidth. 

Obviously,  we  can  apply  the  constant  model  results  for  a  ranging  radar 
directly  to  the  case  of  a  Doppler  radar.  The  ranging  radar  constant  PT 
results  don't  translate  to  Doppler  results  so  easily  because  the  two 
radars  have  opposite  CNR  dependencies  on  Q.  However,  we  can  make  plots 
like  those  in  figures  6.8  and  6.9  for  this  new  CNR  dependence  on  Q. 

Figures  6.10  (a)  and  (b)  plot  the  same  performance  and  contrast  cases  as 
in  figures  6.8  and  6.9  but  for  this  new  Q  dependence.  Now  as  we  decrease 
the  Doppler  velocity  resolution,  the  transmitted  pulse  must  grow  narrower 
and  so  the  IF  bandwidth  must  grow  larger.  This  causes  the  CNR  to  drop. 

The  data  base  used  for  all  the  plots  in  this  chapter  only  goes  down  to 
CNRt  *  10  dB,  so  when  Q  falls  low  enough  that  CNRt  falls  to  below  10  dB, 
we  assume  the  processor  cannot  meet  the  performance  requirements. 
(Realistically,  we  don't  expect  much  out  of  any  processor  at  such  low  CNR 
levels.)  Because  the  CNRfc  values  fall  so  quickly,  figure  6.10  (c)  plots 
the  same  performance  and  contrast  cases  as  in  figures  6.10  (a)  and  (b) 
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Figure  6 . 10 

Doppler  Only-Intensity  Only  Processor  Performance  Trades: 

Constant  P^  Model 

(a)  High  Performance,  High  Contrast 

(b)  Low  Performance,  Low  Contrast 

(c)  High  CNR^0)  Cases 
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except  for  CNR^  *  30  dB  so  we  can  get  more  points  on  the  plot  before 


CNR^  falls  below  10  dB. 


A  drop  in  Q  from  10000  to  1000  reduces  the  system  CNR  by  10  dB.  This 
effect  dominates  the  performance  plots  in  figure  6.10.  This  effect  is 
strong  enough  to  overwhelm  the  much  weaker  underlying  Q  dependences  in  the 
processors  which  we  plotted  in  figures  6.2  (a)  and  6.5  (a).  Both 
processors  exhibit  degraded  performance  for  decreased  Q.  For  any  fixed  N 
value,  however,  the  Doppler  only  processor  always  performs  at  least  as 
well  as  the  intensity  only  processor  in  these  cas^c. 


The  poorer  performance  of  the  intensity  only  processor  at  low  contrasts  is 
still  readily  apparent  in  these  plots.  If  we  have  only  5  dB  of  contrast, 
even  at  CNR^^  =  30  dB  we  need  at  least  N  =  30  to  achieve  the  performance 
standards  with  the  intensity  only  processor.  At  the  same  CNR^0)  value, 
the  Doppler  only  processor  meets  performance  requirements  at  N  =  10  for 
both  performance  and  contrast  cases. 


The  plots  in  figure  6.10  give  us  additional  insight  into  the  effect  of 
changing  bandwidth  over  that  gained  from  the  plots  in  figure  6.9.  In 
figure  6.9,  the  change  in  bandwidth  only  improved  an  already  reasonable 
CNR.  Here,  the  change  (increase)  in  bandwidth  reduces  the  CNR  and  we  can 
see  quite  dramatically  the  effect  the  falling  CNR  has  on  performance. 


We  will  conclude  this  section  and  the  chapter  by  doing  a  brief  analysis  of 
the  range  and  Doppler  case  using  the  thumbtack  ambiguity  function  of 
figure  3.3,  where  T  is  the  duration  of  the  waveform  and  B  is  its 


r>ru»y.w\i-  ~ -r  ’ij.  r a  *>. 


bandwidth.  With  this  case  we  have  a  bit  more  difficulty  being  specific 
about  the  relationships  between  CNR  and  T  and  B  because  we  have  not 
specified  the  waveform  explicitly  and  we  cannot  specify  the  duty  cycle, 
t  /T  in  our  notation,  based  solely  on  the  ambiguity  function.  (Actually, 
the  ambiguity  function  uniquely  determines  the  waveform,  but  there  is  no 
general  procedure  for  inverting  an  ambiguity  function  to  find  the 
corresponding  transmitted  waveform.)  Microwave  radars  often  use  coded 
pulses,  up-chirps,  down-chirps,  and  combinations  of  these  waveforms  [20]. 
We  will  assume  a  duty  cycle  of  1  (or  at  least  a  duty  cycle  independent  of 
T  and  B)  so  the  P^-limited  and  Pp-liraited  cases  are  identical.  This  is  a 
good  assumption  at  least  for  the  CW  chirped  waveforms.  Now  we  will  define 
two  Q  related  values:  Qp  =  FT,  the  number  of  Doppler  bins,  and  =  L^B, 
the  number  of  range  bins,  with  Q  =  QpQ^  (see  figure  3.4).  For  a  unity 
duty  cycle. 
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exactly  as  for  the  ranging  radar  with  constant  Pp.  There  is  another 
consideration,  however,  which  won't  permit  us  to  use  the  ranging  radar 
analyses  directly  when  changing  the  value  of  Q  or  Qp.  When  computing 
performance  numbers,  we  must  use  the  value  Q  =  QpQp  since  it  is  this  total 
Q  which  influences  the  peak  detector  and  Pr(  a ^  ).  However,  when  we 
change  Q  and  examine  the  effect  on  the  performance  through  the  scaled  CNR, 


we  must  scale  only  on  changes  in  Qp  with  no  scaling  for  changes  in  Qp. 


This  brief  analysis  is  meant  only  to  show  how  the  scaling  laws  can  be  used 
in  a  range  and  Doppler  analysis.  For  specific  waveforms,  more  precise 
analyses  are  possible.  However,  we  should  still  expect  to  see  different 
dependences  on  and  Qp  in  the  scaling  laws  while  the  performance 
equations  remain  a  function  only  of  Q. 
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VII.  SUMMARY  AND  CONCLUSIONS 

In  this  thesis  we  have  attacked  the  problem  of  obtaining  analytical 
expressions  for  the  performance  of  laser  radars  as  a  function  of  several 
signal  dimensions  for  multipixel  targets.  To  derive  such  expressions,  we 
proposed  a  reasonable  model  for  real  target?  and  their  environments  in  the 
form  of  a  vertical  target  against  a  sloping  background  at  unknown  range 
and  location.  We  derived  an  expression  for  the  IF  signal  in  the  radar  for 
a  single  pixel  based  on  the  diffuse  speckle  target  model  and  used  the 
components  of  our  target  and  environment  model  to  incorporate  the  range 
and  spatial  information.  We  derived  near  optimal  intensity,  range,  and 
joint  range  and  intensity  measurement  processors  sacrificing  optimality 
for  simplicity  where  necessary  without  over-simplifying  the  model.  We 
were  able  to  derive  performance  measures  for  two  of  the  processors  and 
compute  the  performance  as  a  function  of  the  radar  parameters  and 
resolutions.  Finally,  we  examined  some  specific  cases  of  the  results  and 
explained  physically  what  phenomena  were  at  work  to  give  the  observed 
results. 

What  we  found  qualitatively  was  that  the  range  only  processor  generally 
outperforms  the  intensity  only  processor.  This  comes  as  no  surprise  to 
researchers  in  this  field.  However,  we  were  also  able  to  make 
quantitative  statements,  and  this  is  something  new  to  this  area  of 
research.  For  a  given  set  of  parameters,  we  can  say  which  processor  will 
perform  better  and  by  how  much.  We  can  deduce  the  regions  where  an 
intensity  only  processor  outperforms  a  range  processor  and  by  how  much. 

We  made  quantitative  statements  about  how  sensitive  performance  is  to  any 


of  the  system  parameters  we  investigated  (N,  Q,  CNRt>  and 

Although  we  achieved  our  stated  goal  for  the  thesis,  there  are  several 
areas  where  future  work  could  extend  the  results.  The  most  important  of 
these  areas  involve  relaxing  assumptions  and  making  the  models  even  more 
realistic. 

The  first  assumption  which  we  might  relax  is  the  assumption  the  target 
aligns  exactly  with  a  subframe.  Although  this  assumption  is  probably  the 
most  physically  unappealing  of  all,  it  may  be  the  least  important.  We  saw 
(equation  4.7)  the  system  performance  depends  on  the  performance  for  a 
single  subframe  of  data.  Once  we  have  the  performance  or  statistics  for  a 
single  subframe,  the  overall  system  performance  comes  directly  from 
maximizing  over  the  subframes.  If  we  take  a  window  the  size  and  shape  of 
the  target  and  slide  it  around  the  frame  (like  the  dotted  target  outline 
in  figure  2.5)  collecting  a  subframe  statistic  for  each  possible  window 
location,  we  would  collect  about  MN  subframes  statistics,  each 
statistically  correlated  with  N  other  subframe  statistics.  Consequently, 
the  difference  between  this  model  and  the  target  aligned  with  the  subframe 
model  is  simply  maximizing  over  MN  dependent  random  variables  or  over  M 
independent  random  variables.  Whether  we  maximize  over  independent 
subframes  or  statistically  dependent  subframes  tells  us  nothing  new  about 
laser  radars  per  se.  Furthermore,  the  number  of  subframes  is  independent 
of  any  of  the  radar  parameters  (chapter  6,  section  A).  The  one  factor 
which  does  change  is  the  statistics  of  the  subframes  if  the  target  is  not 
aligned  exactly  with  the  subframes.  However,  in  the  absence  of  noise,  the 
maximum  contrast  in  either  intensity  or  range  space  must  occur  when  the 
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window  aligns  with  the  true  target  location.  All  these  factors 
considered,  it  may  be  more  productive  to  concentrate  on  better 
characterizations  of  the  subframe  statistics  rather  than  better  models  for 
the  subframe  maximization  procedure. 

For  the  subframe  model  we  used,  the  most  important  extension  of  the  work 
is  probably  an  analysis  of  the  joint  inuensity  and  range  processor  in 
chapter  4,  section  F.  This  will  probably  involve  simulation,  but  there 
may  yet  be  further  analytical  steps  an  investigator  could  take. 

The  final  general  area  for  future  work  on  the  models  is  in  the  model  for 
the  target  and  its  environment.  Here  there  are  several  ways  the  models 
could  improve. 

First,  there  are  still  parameters  which  we  should  have  treated  as  unknown 
in  a  completely  general  model.  Two  of  these  are  the  target  size  and  the 
target  shape.  Treating  target  size  as  an  unknown  parameter  will  involve 
adding  N  to  the  set  of  unknown  parameter  in  vector  A  in  chapter  4. 

Assuming  an  unknown  target  shape  will  be  a  bit  more  difficult  to 
incorporate  into  the  generalized  likelihood  ratio  models.  One  potential 
method  is  to  define  a  standard  set  of  shapes  and  maximize  over  the 
discrete  members  of  that  set. 

Second,  if  the  background  is  not  smooth,  our  assumption  of  known  range  to 
background  is  no  longer  valid.  Then  a  ranging  radar  with  very  good  range 
resolution,  lr,  might  make  many  more  errors  because  it  used  the  wrong  bins 
for  counting  the  background  hits,  K.  One  solution  to  this  problem  could 
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be  to  assume  the  background  is  only  reasonably  smooth  and  expand  the  size 
of  the  set  of  background  range  bins.  That  is,  increase  Qg  for  the  same  Q, 
then  perform  a  maximization  over  the  Qg  bins  similar  to  the  target 
maximization  to  get  a  value  for  k.  However,  since  the  correct  bin  number 
changes  from  one  scan  line  to  the  next  and  may  now  also  change  within  a 
scan  line  since  the  background  is  not  smooth,  this  maximization  will  be 
much  more  difficult  to  define  and  perform. 


Third,  both  target  and  background  could  have  non-uniform  reflectivity,  so 
CNRt  and  CNRg  would  be  functions  of  angular  location  (vary  by  pixel).  We 
have  assumed  throughout  the  thesis  the  CNRs  are  independent  of  position: 
the  target  and  background  have  uniform  reflectivity.  Particularly  for  the 
background,  which  is  much  larger  than  the  small  targets  we  used  in  this 
thesis,  this  assumption  is  probably  not  valid  in  any  real  problem  of 
interest. 
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